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FOREWORD 


This  report  is  one  of  a series  of  reports  that  describe  work  performed 
by  Douglas  Aircraft  Company,  McDonnell  Douglas  Corporation,  .3855  Lakewood 
Boulevard,  Long  Beach,  California  90846,  under  the  Windshield  Technology 
Demonstrator  Program.  This  work  was  sponsored :by  the  U.  S.  Air  Force  Flight 
Dynamics  Laboratory,  Wright-Patterson  Air  Force  Base,  under  Contract 
F3361 5-75-C-31 05 , Project  2202/0201  . 

This  report  is  divided  into  three  parts.  Part  1 is  entitled  "Theory 
and  Application",  Part  2 is  entitled  "User's  Manual",  and  Part  3 is  entitled 
"Programming  Manual."  The  principal  investigators  and  authors  were  P.  H.  Denke 
for  Part  1 , G.  R.  Eide  for  Part  2,  and  R.  C.  Morris  for  Part  3. 

Mr.  D.  C.  Chapin,  Capt. , USAF  Ret.,  was  the  Air  Force  Project  Manager 
during  the  conceptual  phase  of  the  work  reported  herein.  Lieutenant  L.  G. 
Moosman  (AFFDL/FEW)  succeeded  Mr.  Chapin  during  the  conduct  of  the  program. 

Mr.  J.  H.  Lawrence  Jr.,  was  the  Program  Director  for  the  Douglas 
Aircraft  Company. 

This  report  was  submitted  to  the  Air  Force  on  7 December  1977,  and 
covers  the  work  performed  during  the  period  July  1975  through  December  1977. 
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SECTION  I 
INTRODUCTION 

In  recent  years  design  verification  of  aircraft  windshield  systems  sub- 
jected to  bird  Impact  has  been  accomplished  primarily  by  test.  A reliable 
analysis  tool  is  needed  to  support  the  design  process  and  reduce  the  amount 
of  testing  required  to  substantiate  the  design.  This  document  gives  detailed 
instruction  in  the  use  of  such  a tool— the  Bird  Impact  Math  Model  Computer 
Program  (IMPACT). 

PROGRAM  FUNCTION 

IMPACT  is  a CDC  oriented  computer  program  for  calculating  the  transient 
dynamic  response  of  a windshield  system  to  loads  resulting  from  bird  impact. 
The  theoretical  development  and  applications,  and  the  program  documentation 
and  implementation  procedure  are  found  in  Parts  1 and  3 (separate  volumes) 
of  this  report.  Appendix  A contains  specific  operational  information  In 
the  form  of  direct  duplication  of  Sections  II  and  III  of  Part  3— Programming 
Manual.  These  sections  are  oriented  to  the  Wright  Patterson  Computing 
Facility  and  should  be  directly  applicable  for  users  at  that  location.  Refer 
to  Part  3— Programming  Manual  for  additional  information  on  implementation 
and  execution  at  other  installations. 

PROGRAM  CAPABILITIES 

IMPACT  solves  the  equations  of  motion  of  a finite  element  representation 
of  the  windshield  system  which  can  Include  multi-layered  transparencies. 

In  fact,  the  analysis  capability  is  not  limited  to  windshield  systems  but 
may  be  applied  to  a finite  element  model  of  any  structure. 

Pre-processing  features  Include  the  capability  to  generate  joint  coordin- 
ates and  element  definition  data  for  a multi-layered,  variable  thickness 
transparency  from  input  data  defining  the  surface  of  the  transparency.  Also 
Included  Is  the  capability  to  generate  model  joint  loads  representing  the 
bird  impact  by  Inputting  the  bird  mass  and  velocity  together  with  model 
geometry  and  Impact  location  data.  In  addition  to  generated  Impact  loads, 
the  user  may  Independently  Input  other  joint  loads  which  may  or  may  not  vary 
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with  time.  Optionally,  the  preprocessing  features  may  be  by-passed,  and 
the  required  card  inout  may  be  provided  by  the  user. 


The  presence  of  initial  thermal  gradients  in  the  structure  may  be 
accounted  for  by  inputting  temperatures  at  the  joint  locations  in  the  finite 
element  model.  Mot  only  are  thermal  deformations  accounted  for,  but  also 
the  material  properties  are  calculated  as  a function  of  temoerature  on  an 
element  by  element  basis.  The  linear  version  of  the  program  accounts 
for  thermal  deformations  in  the  solution.  However,  in  the  current 
version  of  the  nonlinear  program,  the  code  which  incorporates  thermal 
deformations  has  been  deactivated. 

Modal  transformation  of  the  equation  of  motion  is  used  to  reduce  the 
problem  size  to  a practical  limit,  so  that  a stepwise  solution  in  time  may 
be  generated  Involving  a large  number  of  time  increments.  Either  a linear 
or  nonlinear  solution  of  the  incremental  response  may  be  chosen  by  the  user. 

The  linear  solution  is  based  on  linear  small  deflection  theory  and  is  appli- 
cable to  those  problems  which  are  expected  to  respond  in  a linear  fashion. 

Since  It  is  more  economical  than  the  nonlinear  solution,  it  can  be  advan- 
tageous to  use  the  linear  solution  to  perform  comparative  design  studies  of 
a nonlinear  system  followed  by  a final  nonlinear  solution  of  the  chosen 
design. 

The  nonlinear  solution  features  an  Iterative  solution  within  each  time 
Increment  using  the  fictitious  force  approach  to  account  for  geometric 
nonlinearity  with  updating  of  joint  coordinates  and  the  calculation  of  an 
imbalance  correction  to  the  applied  loads  before  proceeding  to  the  next  time 
Increment.  Current  experience  Indicates  that  the  type  and  number  of  modes 
used  In  the  modal  transformation  may  be  very  significant  in  nonlinear  analysis. 
Therefore,  for  accurate  results,  care  and  understanding  must  be  exercised 
In  the  use  of  transformation  modes.  Experience  has  also  shown  that  non- 
linear results  are  very  accurate  for  both  high  and  low  damped  systems  which 
are  small  enough  to  be  analyzed  without  modal  transformation.  The  nonlinear 
solution  Is  designed  to  account  for  geometric  nonlinearity,  Including  large 
deflections,  and  material  nonlinearity.  However,  all  aspects  of  the  material 
nonlinearity  capability  are  not  complete,  so  that  feature  is  currently  not 
available. 
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Both  the  linear  and  nonlinear  solutions  output  a file  of  data  containing 
the  modal  incremental  displacements,  accumulated  displacements,  velocities 
and  accelerations  together  with  element  forces,  stresses  (including  the 
Prandtl-Reuss  equivalent  stresses)  and  strains  for  each  time  point.  This 
file  is  generated  on  magnetic  tape.  Post-processor  capability  allows  the 
user  to  selectively  print  the  modal  response  data,  untransformed  displace- 
ments, joint  coordinates,  and  element  forces,  stresses  and  strains. 

The  linear  or  nonlinear  static  response  of  a system  can  also  be  calcu- 
lated if  the  solution  is  continued  until  the  incremental  response  becomes 
Insignificant.  Static  equilibrium  is  reached  in  the  shortest  time  if  arbi- 
trarily high  damping  is  introduced. 

UNITS 

The  IMPACT  program  has  been  developed  without  dependence  on  units.  The 
user  may  choose  any  system  of  units  he  wishes,  and  he  must  then  be  consistent 
in  the  use  of  those  units  for  all  data  input  to  the  program.  Results  must  be 
interpreted  by  the  same  units  used  for  the  input. 

PROJECTED  IMPROVEMENTS 

The  following  improvements  are  considered  to  be  important  in  the  further 
development  of  IMPACT: 

1)  Improved  Impact  load  generation  Including  pressure  distribution 
within  the  footprint  and  the  effects  of  bird  and  compliant 
target  interaction. 

2)  Automated  mode  extraction  tailored  to  nonlinear  requirements. 

3)  Completion  of  the  material  nonlinearity  capability. 

4)  Addition  of  multi -connect  capability. 

5)  Improved  postprocessor  capability  Including  Interactive  and/or 
batch  mode  plot  capability. 

6)  Improved  user  convenience. 

7)  Addition  of  restart  capability. 
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SECTION  II 


MODELLING  CONSIDERATIONS 

This  section  contains  information  needed  by  the  user  to  create  a 
finite  element  model  of  a windshield  system  which  may  include  a 
multi-layered  transparency  and  support  structure.  It  is  assumed  that 
the  user  is  familiar  with  the  basic  conceDts  involved  in  creating  a 
finite  element  model.  Therefore,  this  discussion  will  be  limited  to 
modelling  procedures  directly  related  to  the  IMPACT  program. 

Instructions  for  preparing  the  data  cards  which  define  the  model  to 
the  program  are  found  in  Section  III  under  INPUT  DATA  REQUIREMENTS,  Data 
Code  Descriptions. 

FINITE  ELEMENTS 

Figure  1 shows  the  finite  elements  which  may  be  used  to  model  a 
structure  for  analysis  by  IMPACT.  Detailed  discussions  and  deriva- 
tions of  these  elements  are  presented  in  "Part  1,  Theory  and  Appli- 
cations", which  is  a separate  volume  of  this  report.  They  are 
called  "lumped  parameter"  finite  elements.  In  the  lumped  oarameter 
concept,  membrane  elements  and  solid  plate  (cell)  elements  are 
assembled  from  bars  which  carry  only  axial  load  and  panels  which 
carry  only  shear.  The  lumping  process  is  automated,  so  user  input 
is  reduced  to  basic  geometric  and  material  data.  The  element  forces 
available  from  the  Postprocessor  are  the  lumped  element  forces  as 
shown  in  Figures  1 and  2,  The  user  may  convert  them  to  stresses  at 
'the  joints  by  multiplying  by  the  appropriate  stress  transform  (see 
Part  1 Appendix  C on  Element  Derivations).  They  can  also  be  resolved 
to  forces  and  moments  at  the  mid-plane  of  a plate  element.  Element 
stresses  and  strains  calculated  by  the  program  and  printed  by  the 
Postprocessor  are  average  values  at  element  mid  points. 


"i  ■■■!"-  

hbcedino  pact  bunk-hot  films? 
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Axial  Bar 


The  bar  element  carries  axial  load  only.  It  is  defined  by  the 
user  on  a data  code  20  card  by  entering  an  identification  number,  two 
joint  numbers  p and  q at  the  ends  of  the  bar,  a material  reference 
number,  and  the  cross-sectional  area.  There  are  two  lumped  element 
forces  as  shown  in  Figure  1.  The  average  stress  and  associated  strain 
are  calculated  at  the  mid  point  of  the  bar. 

Membrane 


The  membrane  element  carries  bi-axial  load  and  shear.  It  is 
defined  by  the  user  on  a data  code  30  card  by  entering  an  identifica- 
tion number,  four  joint  numbers  p,  q,  r,  s at  the  corners  of  the 
membrane  thickness,  and  a stress  orientation  angle.  There  are  nine 
lumped  element  forces  as  shown  in  Figure  1.  The  three  stress  compo- 
nents represent  the  average  condition  of  plane  stress  at  the  mid  point 
of  the  membrane,  and  the  orientation  with  respect  to  the  pc  edge  is 
controlled  by  the  user  input  angle  fpq  (degrees).  In  the  case  of  no 
input  the  default  orientation  will  be  parallel  and  perpendicular  to 
the  pq  edge. 

This  element  may  be  warped  and  may  have  a general  quadrilateral 
shape,  but  accuracy  will  be  best  if  it  is  a parallelogram.  Bi-axial 
Poisson  coupling  is  accounted  for. 

Cell 


The  cell  element  carries  all  types  of  load.  It  is  defined  by  the 
user  on  a data  code  40  card  by  entering  an  identification  number, 
eight  joint  numbers  pQ,  qQ,  rQ,  sQ,  p-j , o-| , r^ , s-j  at  the  corners  of 
the  cell,  a material  reference  number,  and  a stress  orientation  angle, 
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See  Figure  1.  It  is  used  to  represent  a solid  hexahedron  which  should 
be  an  approximate  parallelepiped  for  best  accuracy.  It  is  called  a 
cell  because  the  lumped  parameter  representation  is  a cell  composed  of 
12  bars  and  6 shear  panels.  It  is  basically  a plate  element  where  the 
joints  Pq  through  Sg  are  in  the  upper  surface  of  the  plate  and  the 
joints  p.j  through  s-j  are  in  the  lower  surface.  The  plate  may  be  thick 
or  thin  - joint  coordinates  place  the  joints  in  the  true  plate 
surfaces  thereby  defining  the  thicknesses  at  the  four  corners.  Thus, 
a laminated  transparency  may  be  modelled  realistically  by  a stack  of 
cell  elements  each  having  the  true  thickness  and  properties  of  the 
layer  it  represents. 


If  a string  of  cell  elements  is  used  to  model  a beam  having  a 
cross-sectional  area  of  A and  vertical  and  lateral  bending  moments  of 
inertia  of  Iy  and  1^,  then  the  width  w and  thickness  t of  the  cell 
elements  in  the  model  may  be  calculated  from 


w 


16  \ 1/8 

~) 


t a 


pj3y/8 


(i) 


(2) 


Such  a beam  model  will  have  the  correct  bending  stiffnesses  both 
vertically  and  laterally,  but  axial  stiffness  will  be  sacrificed.  If 
axial  stiffness  is  more  Important  than  lateral  bending,  the  following 
equations  may  be  used  to  satisfy  vertical  bending  and  axial 
stiffnesses. 
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The  six  stress  components  shown  in  Figure  1 represent  the  average 
state  of  stress  at  the  mid-points  of  the  uoper  and  lower  surfaces. 
Therefore,  twelve  average  stresses  and  strains  are  calculated  for  a 
cell  element.  The  x and  y comoonents  are  in  the  plane  of  the  surface, 
and  the  x component  is  oriented  with  respect  to  the  PqPq  edge  by  the 
user  input  angle  fpg  (degrees).  Figure  2 shows  the  lumped  element 
forces  for  a cell  element.  The  lumped  elements  are  located  at  a 
reduced  depth  of  t/v^  in  order  to  oroduce  the  correct  bending 
stiffness,  Tri -axial  Poisson  coupling  is  accounted  for. 

Point  Mass 

A fourth  element  not  shown  in  Figure  1 is  described  here.  The 
point  mass  element  represents  a local  mass  attached  to  the  model  at 
a specified  joint.  It  is  defined  by  the  user  on  a data  code  50  card 
by  entering  an  Identification  number,  the  joint  number  where  the  mass 
Is  located,  and  its  contributions  in  the  global  X,Y,Z  degrees  of  free- 
dom. These  contributions  are  added  to  the  appropriate  diagonal  ele- 
ments of  the  structural  mass  matrix.  This  element  is  provided  for 
the  purpose  of  Including  non -structural  mass  In  the  finite  element 
model . 
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GRID  DETAIL  AND  GEOMETRY 


Model  development  usually  begins  with  a drawing  or  sketch  of  the 
structure  on  which  the  user  lays  out  a grid  representing  subdivision 
of  the  structure  into  finite  elements.  Surface  qrids  and  cross- 
sectional  grids  define  the  three  dimensional  model.  Scale  drawings 
are  helpful  for  visualizing  finite  element  shapes  and  for  designing 
the  grid  so  that  important  structural  features  are  retained  while 
varying  the  grid  spacing  so  that  model  detail  is  provided  in  critical 
regions  of  the  structure.  Intersections  of  grid  lines  are  called 
joints,  and  model  geometry  is  defined  in  terms  of  joint  coordinates. 


Generation  of  joint  coordinates  is  probably  the  most  difficult  and 
time  consuming  task  in  creating  a finite  element  model  of  a complex 
structure.  Some  joints  may  be  located  exactly  at  control  points  for 
which  coordinates  are  available  on  a drawing.  Sometimes  joint 
coordinates  may  be  scaled  from  drawings.  Coordinates  for  joints  in 
the  surface  of  a structure  are  usually  calculated  from  loft  data  while 
interior  joints  are  located  by  offset  from  the  lofted  surface.  Some 
coordinates  have  to  be  calculated  from  equations  that  satisfy  section 
property  requirements  in  local  structural  elements.  Some  instal- 
lations may  have  interactive  graphics  or  other  computer  aided 
techniques  to  help  the  user  in  this  process. 

Transparency 

In  a windshield  system  the  transparency  may  be  a single  layer  or 
a multi-layered  laminate.  In  either  case  the  user  has  to  lay  out  the 
grid  and  calculate  the  joint  coordinates  for  the  outer  surface  only. 

He  may  then  use  the  Laminate  Generator  program  to  generate  interior 
joint  coordinates.  See  Section  III  under  INPUT  PREPROCESSING,  Laminate 
Generator.  The  grid  mesh  should  be  fine  enough  so  that  the  dynamic  response 
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of  the  structure  can  be  adequately  represented.  Because  of  model  size 
limitations  (1200  joints)  and  if  there  are  many  layers  in  a laminated 
transparency,  it  is  sometimes  difficult  to  provide  the  desired  degree 
of  detail.  Compromise  may  be  necessary  in  the  form  of  coarser  grid 
in  areas  remote  from  the  imoact  point,  although  careful  consideration 
should  be  given  as  to  whether  critical  response  might  occur  at 
locations  other  than  the  impact  point.  Model  size  can  also  be 
decreased  by  taking  advantage  of  structural  symmetry  and  by  modelling 
less  of  the  support  structure.  The  transparency  will  have  some  joints 
in  common  with  the  supDort  structure.  Therefore,  planning  of  the 
transparency  grid  must  be  coordinated  with  important  features  of  the 
support  structure. 

Support  Structure 

As  stated  above,  the  grid  layout  for  the  support  structure  is 
established  in  coordination  with  that  of  the  transparency.  The  extent 
and  degree  of  detail  to  be  included  in  the  model  depends  on  consid- 
erations of  factors  such  as  model  size  limitations  and  location  of 
impact  point.  Impact  on  or  near  the  transparency  support  members 
would  reouire  more  modelling  detail  to  be  concentrated  in  the  support 
members  than  impact  in  the  center  of  the  transparency. 

Extruded,  formed  or  built-up  framing  members  around  the  trans- 
parency are  usually  modelled  as  beams  composed  of  cell  elements.  Load 
paths  connecting  the  framing  members  to  surrounding  shell  structure 
and  a portion  of  the  shell  structure  itself  may  be  realistically 
represented  by  bar,  membrane  and  cell  elements.  Joint  coordinates 
are  calculated  from  loft  data  and  from  drawing  information. 
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Loads 


Loads  are  applied  to  the  model  at  the  joints  in  the  global  X, 

Y,  Z degrees  of  freedom.  The  calculation  of  bird  impact  loads  at 
specified  joints  is  described  under  "Loads  Generator".  They  are  intro- 
duced in  the  input  stream  and  converted  to  transformed  loads  as 
described  under  "MODAL  TRANSFORMATION". 

ATTACHMENTS 

In  the  overall  model,  the  attachments  between  the  transparency  and 
the  support  structure  must  be  simplified.  A single  shear  bolted 
attachment  may  be  represented  as  a simply  supported  or  hinged  connec- 
tion along  a single  line  of  cotnnon  joints  such  as  joints  16,  39  and  62 
in  Figure  11.  A double  shear  attachment  could  be  represented  by  adding 
membrane  elements  to  connect  edge  13,  36,  59  to  edge  22,  45,  68.  In  a 
laminate  having  structural  plies  separated  by  softer  interlayer  mate- 
rial the  presence  of  bushings  and  bolts  through  the  layers  can  be 
represented  by  specifying  stiff er  shear  properties  for  the  interlayer 
elements  along  the  edges  of  the  laminate.  The  user  can  approximate 
the  stiffer  properties  by  hand  analysis. 

If  necessary  the  user  can  perform  a more  detailed  study  of  laminate- 
bolt-support  structure  Interactions  by  making  a finite  element  model 
of  the  local  attachment  region.  Important  features  such  as  bolt,  bush- 
ing, and  portions  of  the  laminate  and  support  structure  can  then  be 
Individually  modelled  in  much  more  detail  than  in  the  overall  model  of 
the  system.  Response  of  this  localized  attachment  model  to  typical 
impact  loads  can  be  studied  to  improve  understanding  of  the  interactions 
of  the  component  parts  and  to  gain  insights  in  how  to  represent  them 
realistically  in  the  larger  model. 
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CONSTRAINTS 


A joint  may  be  constrained  (fixed)  against  displacement  in  one, 
two  or  three  specified  directions.  Constraint  in  one  direction  still 
allows  the  joint  to  be  displaced  in  a plane  normal  to  the  specified 
direction.  Constraint  in  two  directions  fixes  the  joint  against 
movement  in  the  plane  containing  ttie  two  specified  direction  vectors 
but  still  allows  the  joint  to  be  displaced  normal  to  the  plane.  Con- 
straint in  three  directions  fixes  the  joint  in  space  and  allows  no 
displacement  in  any  direction.  Constraint  data  is  input  on  data 
code  5 cards. 

Constraints  are  used  to  support  the  finite  element  model  in  a 
stable  manner  and  to  approximate  boundary  conditions  at  the  edges  of 
the  model  where  structural  continuity  needs  to  be  represented.  For 
example,  on  a syrnnetrlcally  loaded  structure  the  correct  boundary  con- 
ditions at  an  edge  which  lies  in  a plane  of  structural  symmetry  are 
provided  by  constraining  all  joints  in  the  plane  of  symmetry  in  the 
degree  of  freedom  normal  to  the  plane.  Constraints  decrease  the 
problem  size,  since  each  constraint  eliminates  one  degree  of  freedom. 

MODEL  DIAGRAMS 


The  user  should  prepare  diagrams  of  the  finite  element  model  showing 
grid  lay-outs  and  numbering  sequences  of  joints  and  elements.  If  a 
computer  aided  drawing  or  plotting  capability  is  available  and  can  be 
Interfaced  to  read  model  definition  data  (or  even  to  help  create  It), 
then  scale  drawings  may  be  produced  showing  the  model  from  advantageous 
angles.  Often  non-scale  diagrams  are  also  used  to  clarify  model 
details.  Figures  3 through  8 are  example  diagrams  of  a large  model 
having  910  joints,  28  bars,  143  membranes,  460  cells  and  4999  degrees 
of  freedom. 
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Aft  Frame 


Figure 


Joint  Numbering 

The  order  in  which  joints  are  numbered  is  important.  Decompo- 
sition of  the  stiffness  matrix  is  based  on  a wave  front  technique  which 
is  most  efficient  when  the  wave  front  is  small.  If  the  wave  front  is  too 
large  it  may  be  impossible  to  complete  the  decomposition  because  of  pro- 
hibitive computer  core  requirements.  The  size  of  the  wave  front  is  related 
to  the  joint  numbering  sequence. 

The  wave  front  can  be  minimized  by  numbering  the  joints  in  a 
systematic  and  continuous  manner,  section  by  section,  from  one  end  of 
the  model  to  the  other.  For  this  purpose  a section  may  be  visualized 
as  a slice  through  the  model  completely  severing  it  (except  for  the 
first  and  last  slices),  and  selected  in  such  a way  that  the  number  of 
joints  in  the  section  is  minimized.  Section  Interfaces  need  not  be 
planes.  The  number  of  joints  in  a section  may  vary  from  one  section 
to  the  next.  The  wave  front  size  will  vary  accordingly  during  the 
decomposition  and  will  be  approximately  equal  to  6 times  the  number 
of  joints  in  a section.  The  maximum  wave  front  determines  the  core 
requirements.  Reference  1,  Volume  V Supplement  II,  gives  a thorough 
treatment  of  wave  front  processing  and  optimum  joint  numbering. 

Consideration  should  also  be  given  to  the  numbering  order  within 
a section.  If  the  model  includes  a multi-layered  transparency  it  Is 
usually  best  to  number  the  joints  in  a stack  sequentially  from  outer 
to  inner  surface.  See  "Laminate  Generator"  for  definitions  of  joint 
stacks  and  joint  numbering  increment.  Figure  11  shows  a typical  number- 
ing scheme  where  the  first  section  contains  joints  1 through  23,  the 
second  contains  joints  24  through  46,  and  the  third  contains  joints 
47  through  69.  In  the  model  diagrams  of  Figures  4 through  8 a trans- 
verse section  through  the  three  layered  transparency  and  side  frame 
exhibits  a similar  numbering  sequence.  In  the  vicinity  of  the  forward 
frame  (Figure  4)  the  scheme  has  been  varied  somewhat,  but  in  general 
the  sections  are  transverse  cross  sections  through  the  model . 
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Element  Numbering 

Each  of  the  four  types  of  finite  elements  (bar,  membrane,  cell, 
and  point  mass)  has  its  individual  numbering  sequence.  In  Figures  4 
through  8 bars  are  numbered  1 through  28,  membranes  are  1 through  143, 
and  cells  are  1 through  460.  There  are  no  special  rules  governing 
the  numbering  of  bar,  membrane  and  point  mass  elements.  Any  system- 
atic sequence  that  the  user  considers  convenient  may  be  used.  In  num- 
bering cell  elements  the  user  must  consider  Laminate  Generator  operations 
and  number  the  transparency  surface  elements  so  that  inner  layer  element 
numbers  may  be  derived  by  adding  a constant  numbering  increment.  Refer 
to  the  definitions  and  discussions  under  "Laminate  Generator"  for 
further  information.  In  Figure  8 the  surface  cell  elements  are 
numbered  1 through  133,  the  second  layer  134  through  266,  and  the 
third  layer  267  through  399.  Therefore,  the  numbering  increment  from 
one  layer  to  the  next  is  133. 

MATERIAL  PROPERTIES 

Table  5 lists  the  ten  properties  which  must  be  input  for  each  material 
used  In  the  model.  Reference  3 contains  the  results  of  specialized  testing 
to  determine  the  mechanical  properties  for  polycarbonate  materials.  Standard 
handbooks  are  used  for  metals. 

Ramberq-Osgood 

The  IMPACT  program  uses  a Ramberg-Osgood  representation  of  the  true 
stress/true  strain  curve,  see  Figure  9.  The  analytic  expression  of  the 
curve  is: 


The  program  uses  input  values  for  E,  EA,  aA,  CTp  and  ep  to  calculate  the 
parameters  n,  F and  3^,  These  parameters  are  required  for  material  non- 
linearity considerations.  As  mentioned  In  the  INTRODUCTION,  the  material 
i nonlinearity  code  Is  not  complete,  but  these  basic  calculations  are  coded. 
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Rupture  point 


E 3 slope  tangent  to  Ramberg-Osgood  curve  at  o’  3 0,  7 3 0. 

E 3 slope  of  linear  part  of  curve  (Young's  modulus). 

Ea  * secant  modulus  at  point  A,  approximate  yield  point. 

aL,  7l  3 true  stress  and  strain  at  elastic  limit. 

aA,  7a  3 true  stress  and  strain  at  approximate  yield  point, 

ar.  eV.  = true  stress  and  strain  at  rupture. 

P r 

Figure  9.  Ramberg-Osqood  True  Stress/True  Strain  Curve. 
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and  the  program  expects  a complete  set  of  material  properties  for  each 

material.  Reference  3 provides  properties  1 through  9 for  polycarbonate 

and  interlayer  materials.  For  metals,  assume  that  point  A in  Figure  9 

is  the  intersection  of  the  standard  .002  offset  line  with  the  stress/strain 

curve  at  the  yield  stress  F.  . Then 

**y 


<7fl  On 

-A  • -jr  + .002 


Combi 


ining  equation  6 and  7 and  solving  for  E^  gives 


E.=  (— ^ ) E 
A \ .002E  + F.  I 


From  Reference  3 the  expressions  for  true  stress  and  true  strain  are  evaluated 
at  the  rupture  stress  F^  and  the  rupture  strain  e to  give: 


1 + e 


4) 


€r  = in  (1  + e) 


Property  number  10,  the  damping  to  stiffness  ratio  will  probably  have 
to  be  estimated  for  metals,  unless  test  data  is  available.  For  polycarbonate 
and  Interlayer  materials  values  for  the  damping  ratio  are  available  from 
Reference  4,  where  a method  based  on  combined  testing  and  analytical  pro- 
cedures was  used  to  determine  damping  ratios  for  several  materials. 

It  Is  sometimes  desirable  to  specify  elastic  material  properties.  For 

example,  glass  and  some  interlayer  materials  exhibit  negligible  plastic 

characteristics.  In  such  a case  the  user  may  establish  the  values  for  E 

and  5 and  set 
r 

Ea  ■ E,  5.  ■ = 4 01) 
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The  program  will  recognize  this  as  a special  situation,  and  as  a result 
some  of  the  calculations  oriented  toward  material  nonlinearity  (plasticity) 
will  be  eliminated.  The  user  may  recognize  that  since  material  non- 
linearity capability  is  not  fully  operational  anyway,  he  may  use  this 
elastic  option  as  a way  to  simplify  input  until  such  time  as  material 
nonlinearity  is  completed. 

Strain  Rate  Dependence 

The  program  does  not  automatically  account  for  continually  changing 
strain  rates  in  all  elements  of  the  model.  This  is  thought  to  be  pro- 
hibitive. Instead  the  user  estimates  the  expected  strain  rate  and  inputs 
material  properties  that  are  consistent  with  the  estimated  strain  rate. 

Temperature  Dependence 

If  values  of  the  material  properties  are  available  at  more  than  one 
temperature  the  user  will  have  an  indication  of  temperature  dependence. 

As  explained  in  Section  III  under  INPUT  DATA  REQUIREMENTS,  Data  Code 
Descriptions,  Data  Code  10,  as  many  temperature  and  property  value  pairs 
as  are  available  may  be  input  for  each  proDerty.  The  program  then  uses 
the  input  temperature/value  pairs  in  a least  squares  curve  fitting  pro- 
cedure to  calculate  a specified  number  of  coefficients  for  a polynomial 
equation  which  defines  the  property  as  a function  of  temperature.  An 
alternate  option  allows  the  user  to  input  the  coefficients  of  the  poly- 
nomial equation  directly  if  he  so  desires.  The  first  coefficient  is  the 
value  of  the  property  when  there  is  no  temperature  dependence.  Subsequent 
coefficients  introduce  temperature  effects  In  terms  of  powers  of  the  tem- 
perature of  the  finite  element  being  processed.  The  program  calculates 
finite  element  temperatures  from  joint  temperature  data  which  the  user 
inputs  on  data  code  2 or  3 cards.  Note  that  when  no  temoerature  dependence 
is  being  specified  the  second  option  should  be  used  to  input  only  the  first 
coefficient  which  is  then  equal  to  the  value  of  the  property. 


LIMITATIONS  AND  COMPUTING  TIME 


Model  size  limitations  are  given  in  Table  1.  In  general  these 
limits  should  allow  the  user  to  model  a structure  with  as  much  detail 
as  practicable  from  the  standpoint  of  computing  time. 


TABLE  1.  MODEL  SIZE  LIMITATIONS 


ITEM 

MAXIMUM 

NUMBER 

Joints 

1200 

Direction  Cosines 

50 

Constraints  * 

600 

Materials 

20 

Element  Forces 

20000 

Edge  Degrees  of  Freedom 

3000 

In  the  discussion  under  "Joint  Numbering"  it  was  pointed  out  that 
the  wave  front  during  decomposition  of  the  stiffness  matrix  should  be 
kept  as  small  as  possible  to  minimize  computer  core  requirements. 
Table  2 gives  approximate  correspondence  between  wave  front  and  core 
requirement. 


TABLE  2.  WAVE  FRONT  AND  COMPUTER  CORE 


WAVE  FRONT 

CORE  REQUIREMENT 
(CM) 

140 

100 ,000 

280 

200 ,000 

370 

300,000 

Table  3 is  presented  to  give  the  user  a basis  for  estimating  com- 
puter time.  The  tabulated  values  are  actual  run  times  on  the  CDC  Computer. 
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TABLE  3.  CDC  COMPUTER  TIME,  SECONDS 


[ 


Joints 

910 

599 

16 

Bars 

28 

0 

2 

Model 

Membranes 

143 

0 

2 

Size 

Cells 

460 

215 

2 

Deg.  of  Freedom 

4999 

3042 

76 

Modes 

30 

30 

10 

Wave  Front 

253 

172 

— 

Laminate  Generator 

30  CPU 

11  CPU 

Loads  Generator 

15  CPU 

14  CPU 

Small  model. 

Initial 

Generator 

537  CPU 

1098  I/O  . 

245  CPU 

401  I/O 

all  steps 

combined 

in  one 

computer 

run  for  10 

Decomposition 

4066  CPU 

7324  I/O 

1000  CPU 

1800  I/O 

time  incre- 

Modal Transformation 

5240  CPU 

8320  I/O 

2000  CPU 

3200  I/O 

ments. 

30  CPU 

Data  Files  for  Inc  Sol 

600  CPU 

1500  I/O 

250  CPU 

600  I/O 

80  I/O 

Linear 

Incremental  Sol 

750+20/ Inc 
1250+35/Inc 

292+9/Inc 

500+15/Inc 

Nonlinear  Incremental  Sol 

— 

100+430/ Inc 
100+280/ Inc 

Postprocessor 

37  CPU 

474  I/O 

H 

KEY -PUNCH  DATA  FORMS 


The  following  seven  pages  are  suggested  forms  on  which  the  user  may 
list  model  definition  data  in  order  to  have  it  keypunched  to  data  cards. 

These  forms  are  consistent  with  the  descriptions  given  in  SECTION  III  under 

"Data  Code  Descriptions."  Note  that  the  same  data  card  format  may  be 
applicable  to  more  than  one  data  code.  For  example,  data  codes  1 through  5 

all  require  the  same  data  card  format.  The  user  may  copy  these  forms  as 

needed,  or  he  may  use  them  as  a guide  to  design  his  own  forms. 


RIGHT  ADJUSTED  INTEGER  DATA | FLOATING  POINT  DATA 


ELEMENT  DEFINITION  DATA 

RIGHT  ADJUSTED  INTEGER  DATA  FfLOATING  POINT  DATA  TYPE  OF  DATA 


I in  ii  rrrrr  mr  rr  t ti  i r nr 


SECTION  III 
USING  THE  PROGRAM 


The  use  of  the  IMPACT  computer  program  is  described  in  this  section. 
FUNCTIONAL  FLOW  CHART 

Figure  10  shows  the  sequential  steps  the  user  follows  In  the  application 
of  the  IMPACT  program.  User  generated  card  Input  data  Is  primarily  finite 
element  model  definition  data  but  also  Includes  all  the  rest  of  the  card 
Input  that  will  be  required  during  the  course  of  the  analysis.  The  prepara- 
tion of  input  data  cards  is  described  in  detail  in  the  subsection  INPUT 
DATA  REQUIREMENTS. 

The  CDC  UPDATE  procedure  may  be  used  to  input  the  model  definition  card 
deck  and  store  It  in  an  UPDATE  data  file.  UPDATE  may  be  used  also  to  edit 
a model  data  file  and  to  input  model  data  to  the  Laminate  Generator,  Loads 
Generator,  or  Initial  Generator  programs.  Optionally,  the  user  may  input 
the  model  definition  cards  directly  to  these  programs. 

The  Laminate  Generator  program  is  designed  to  read  a partial  model 
definition  deck  and  generate  the  missing  data  to  form  a complete  model 
definition.  The  Loads  Generator  program  extracts  the  necessary  information 
from  a model  definition  deck  to  generate  bird  impact  loads.  Functions 
identified  to  this  point  are  referred  to  as  "preprocessing"  functions,  and 
the  modules  involved  are  collectively  termed  the  Preprocessor. 

The  Initial  Generator  reads  a complete  model  definition  deck  and  generates 
the  structural  matrices  and  information  lists  required  in  subsequent  analysis 
steps.  Structural  mass  generation,  stiffness  decomposition,  vibration  mode 
extraction  and  modal  transformation  are  performed  by  the  FORMAT  program, 
which  is  described  In  Reference  1.  FORMAT  is  also  used  to  produce  properly 
ordered  data  files  for  Input  to  the  Incremental  solution.  The  Incremental 
Solution  code  calculates  the  modal  response  and  element  forces,  stresses  and 
strains  for  the  specified  time  points,  and  finally  the  post  processor  prints 
out  selected  Information  for  use  in  evaluating  results. 
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Figure  10.  IMPACT  Functional  Flow  Chart. 


44 


Detailed  discussions  of  the  functions  shown  in  Figure  ' / are  pre 
sented  In  the  following  subsections.  Appendix  A contains  an  example 
deck  setup  including  JCL. 


INPUT  DATA  REQUIREMENTS 

User  generated  card  input  data  is  described  in  this  subsection. 


Data  Code  Descriptions 

Model  definition  data  for  the  IMPACT  program  Is  Identified  by  a 
data  code  entered  In  card  columns  1 and  2 of  each  card.  The  data 
code  is  an  integer  value  In  the  range  1 through  99.  The  following 
pages  give  the  data  code  designation,  type  of  data,  input  card  format, 
and  descriptions  of  the  data  to  be  entered  in  each  field  of  the  input 
card.  Card  diagrams  for  each  data  code  show  general  headings  which 
Indicate  whether  the  data  to  be  entered  Is  Integer  or  floating  point, 
then  the  card  fields  are  marked  off  with  a pre-entered  data  code  In 
card  columns  1,  2 and  field  heading  symbols  shown  In  applicable  fields 
(unused  fields  are  crossed  out),  and  finally  the  beginning  card  column 
is  Identified  at  the  left  side  of  each  field.  In  some  cases,  the 
ending  card  column  for  a field  is  also  shown  for  clarity.  Definitions 
of  the  field  heading  symbols  are  listed  following  the  card  diagram. 

A data  code  deck  consists  of  as  many  cards  as  necessary  to  Input 
the  type  of  data  Identified  by  a given  data  code  followed  by  an  end 
card  having  the  data  code  and  four  nines  entered  In  card  columns  1 
through  6.  A model  definition  deck  Is  made  up  of  applicable  data 
code  decks  stacked  In  the  order  of  Increasing  data  codes.  A complete 
model  definition  deck  as  required  by  the  Initial  Generator  program 
must  Include  data  codes  decks  Indicated  In  Table  4.  A data  code  deck 
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may  be  nothing  more  than  the  end  card  If  that  type  of  data  is  not 
needed  to  define  a particular  model,  but  the  end  card  must  be  present 
as  a minimum.  Note  that  some  data  code  decks  are  optional.  However, 
a complete  model  definition  deck  will  not  be  output  from  the  Laminate 
Generator  unless  all  the  optional  data  code  decks  are  Input. 


TABLE  4.  MODEL  DEFINITION  DATA  CODE  SUMMARY 


OATX 

COOE 

TYPE  OF  DATA 

— 1 lil  M il  1 ■ 1 T i l1  ini  1 

1 

Constants 

OPT 

OPT 

REQ 

2 

Joint  Coordinates 

REQ,  GEN 

REQ 

3 

Joint  Temperatures 

OPT 

OPT 

4 

Direction  Number  Definitions 

OPT 

t 

f 

5 

Joint  Constraint  Definitions 

OPT 

REQ 

6 

Surface  Normal  Vectors 

GEN 

OPT 

7 

Laminate  Definition  Data 

REQ 

\ 

- 

8 

Laminate  Variable  Thickness  Data 

OPT 

OPT 

9 

- Not  Used  - 

- 

■ 

■B 

10 

Material  Properties 

OPT 

OPT 

REQ 

11 

Bird  Impact  Definition  Data 

OPT 

REQ 

OPT 

• 12 

Impact  Footprint  Joint  List 

OPT 

REQ 

OPT 

13-19 

- Not  Used  - 

« 

• 

■ 

20 

Axial  Bar  Element  Definition 

OPT 

OPT 

REQ 

21-29 

- Not  Used'  - 

- 

- 

30 

Membrane  Element  Definition 

OPT 

OPT 

REQ 

31-39 

- Not  Used  - 

- 

• 

■» 

40 

Cell  Element  Definition 

REQ,  GEN 

OPT 

REQ 

41-49 

- Not  Used  - 

• 

« 

- 

SO 

Point  Mass  Element  Definition 

OPT 

OPT 

REQ 

51-99 

- Not  Used  - 

- 

! 

REQ  ■ required 
OPT  ■ optional 
GEN  ■ generated  by  program 

REQ,  GEN  ■ basic  data  required  to  be  Input,  additional  data 
Is  generated  by  program 

- ■ must  not  be  present 


L 
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Data  Code  1,  Constants 

User  inputs  two  data  cards  and  one  end  card. 


i 


I 


Card  1 

a » base  temperature  Tg  (such  as  room  temperature)  for  the 
analysis 

b * mechanical  equivalent  of  heat  J (example:  J = 4.184  x 107 
ergs/calorie) 

Card  2 

a = 1.0  x 106 

= bar  compliance  coefficient  (part  of  material  nonlinearity 
capability  for  which  code  is  incomplete,  but  a value  for 
"a"  is  required) 

b ■ 1.0  x 10-12 

= stiffness  suppression  coefficient  (eliminates  insignificant 
values  in  stiffness  matrix) 

Note:  Values  for  "a"  and  "b"  are  currently  required  to  be  input  by 
user,  but  eventually  should  be  built  into  code. 


i 
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Data  Code  2,  Joint  Coordinates 


User  inputs  one  card  for  each  joint  and  one  end  card.  See  pages  11  and 
76  for  discussions  of  joint  coordinates. 


RT  ADJ  INTEGER 

FLOATING  POINT 

2 

JT 

* 

1 

| 

X 

Y 

Z 

TJ 

I 

m 

m 

19 

34 

49 

mm 

1 

SB 

! E 

rn 

nd  C 

ard 

1 

JT  ■ joint  identification  number  ( <9999) 

X,  Y,  Z * X,  Y,  Z coordinates  of  the  joint  In  the  right  hand 
cartesian  coordinate  system  of  the  model 
Tj  ■ temperature  of  the  structure  at  the  joint  location.  (This 
temperature  implies  thermal  deformations  if  different  from 
the  base  temperature  specified  in  data  code  1.) 
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Data  Code  3,  Joint  Temperature 

This  data  overrides  temperatures  previously  Input  on  data  code  2 cards. 
The  user  Inputs  data  cards  as  necessary  and  one  end  card. 


JTa  ■ joint  number 

JTb  » optional  entry  of  joint  number  to  define  a range  of  joint 
nunfcers  JTa  to  JTb,  Inclusive 

Tj  » temperature  of  the  structure  at  joint  JTa  or  at  all  joints 
In  the  range  JTa  through  JTb,  If  JTb  Is  Input 
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Data  Code  4,  Direction  Number  Definitions 

User  Inputs  one  card  for  each  direction  number  and  one  end  card. 


RT  ADJ  INTEGER 

FIOATIN'G  POINT 

1 

JTn 

X 

*Y 

h 

m 

3 

15 

mmmm 

34 

49 

64  72! 

n 

0 

p999 

En 

i Ca 

rd 

I 

DIR  ■ direction  Identification  number  (1,2,3,...) 

Option  1:  Enter  joint  numbers  JTm  and  JTn  to  define  the  direction 
parallel  to  the  vector  JTm  to  JTn 

Option  2:  Enter  direction  vector  components  i x,  £y,  i ^ (need  not  be 
normalized)  to  define  the  direction 

Note:  This  data  provides  a table  of  direction  numbers  which  can  be 
referenced  to  define  the  direction  of  joint  constraints  in  data  code 
5.  It  is  particularly  convenient  when  a large  number  of  constraints 
are  oriented  in  the  same  direction. 
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Data  Code  5,  Joint  Constraint  Definitions 


The  joint  will  be  constrained  against  displacement  in  the  given  direction. 
User  inputs  one  card  for  each  constraint  and  one  end  card.  See  page  14 
for  discussion  of  constraints. 


RT  ADJ  INTEGER 

FLOATING  POINT 

5 

CON 

JTm 

JTn 

DIR 

iy 

m 

3 

7 

11 

15 

19 

34 

49 

64  72 

l 

11 

E 

id  C 

ard 

CON  = constraint  identification  number  (1,2,3,...) 

JTm  * number  of  the  joint  to  be  constrained 

Enter  joint  number  JTn  to  define  the  constraint 
direction  parallel  to  the  vector  JTm  to  JTn 

Enter  a direction  number  DIR  to  reference  a direction 
defined  in  Code  4 data 

Enter  direction  vector  components  i y,  (need  not 

be  normalized)  to  define  the  constraint  direction 


Option  1: 


Option  2: 


Option  3: 


i 


L 
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Oata  Code  6,  Surface  Normal  Vectors,  reference  Figure  11. 

This  data  Is  generated  by  the  laminate  generator  program  but  may 
be  edited  or  Input  by  cards. 


rr 


RT  ADJ  INTEGER 

FLOATING  POINT 

JTs 

| 

nx 

nY 

nz 

As 

N 

3 

m 

19 

34 

49 

£4  72 

■ 

I 

pi 

i 

nd  C 

i.rd 

JTs  ■ joint  number  of  a joint  on  the  surface  of  the  laminate 

nX*  ’Y’  nZ  ■ X,  Y,  Z components  of  a unit  vector  normal  to  the 
surface  at  joint  JTs  and  pointing  Into  the  laminate 

As  ■ surface  area  of  a section  of  the  laminate  associated  with 
joint  JTs,  reference  Figure  11 


Data  Code  7,  Laminate  Definition  Data 

User  inputs  cards  as  necessary  to  define  a multi-layered  laminate 
followed  by  one  L+99  card  and  one  end  card.  (See  pages  75  through 
80  for  discussion  of  laminate  data. 


blank  for  2 s 0 and  t.  = t 

i o 


7^999|  End  Card 

L * laminate  identification  nurrtoer  (any  multiple  of  100  <9900) 

2 * layer  number  (0,  1,  2,  3,  ••.) 

M i * material  identification  number  for  layer  2 (M^  is  a reference 
to  code  10  data) 

■ joint  numbering  Increment  for  layer  i (J^  = JQut^  + ) , see  page  76 

A * element  numbering  increment  for  layer  2 (E^  * E^_1  + AE^)  , see  page  78 

t|  * thickness  of  layer  2 

» surface  offset  t„  when  2*0 
o 


Note:  A model  of  a windshield  system  could  contain  more  than  one  lamin- 
ated transparency  (windshield,  side  window,  etc.).  The  laminate  identifi- 
cation number  L provides  for  identifying  each  laminate  by  a user  assigned 
number  as  defined  above. 
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Data  Code  8,  Laminate  Variable  Thickness  Data 

User  Inputs  cards  as  necessary  to  define  variable  thickness  layers  In  a 
multi-layered  laminate  followed  by  one  end  card.  See  page  79. 


[2E2MSE 

FLOATING  POINT 

n 

1 

JTa 

JTb 

*1 

*2 

(3 

‘4 

*5 

■ 

if 

3 

■ 

21 

31 

41 

51 

61  7q 

M 

1 

1 

999S 

End  Card 

JTa  ■ surface  joint  number  below  which  layer  thicknesses  t],t2, 
t3,  are  to  be  used  to  override  layer  thicknesses 
specified  In  Code  7 data 

JTb  ■ surface  joint  number  which.  If  present,  defines  a range  of 
surface  joint  numbers  JTa  through  JTb  below  which  the  layer 
thicknesses  will  be  effective 

Note:  If  there  are  more  than  six  layers  In  the  laminate,  enter 
tj  through  on  a second  card  with  card  columns  1 through  10 
duplicated  from  the  first  card.  Continue  this  process  as  necessary. 
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Data  Code  10,  Material  Properties 

A set  of  cards  M through  M+99  is  input  for  each  material  used  in  the 
structure.  One  end  card  follows  the  last  set.  See  page  29  for  discussion. 


ALPHAMERIC 


Material  name  and  descriptive  information 


FLOATING  POINT 


Card  col  72  non-blank  for  continuation 
No  restrictions  on  symbol  choice. 


109999  End  Card 


M * user  assigned  material  identification  number  (100,  200,  ...) 
which  can  be  referenced  in  element  definition  cards  to  specify 
the  element  material 

P ■ property  Identification  number  (1,  2,  3,  ...) 

Either  of  two  options  may  be  used  to  Input  each  material  property  as 

a function  of  temperature.  Both  options  result  in  coefficients  a.. 

of  a polynominal  equation  for  the  value  of  the  property  VD  as  a 

2 r 

function  of  temperature  T (Vp  * aQ  + a^T  + a£T  + ...) 

Option  1 

NC  * number  of  polynominal  coefficients  to  be  computed  by  the  program 
from  the  given  (T^ , V^)  pairs 

< number  of  (T^,  V^)  pairs  input 
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Data  Code  10  (Cont.) 


T^,  Vi  ■ pairs  of  values  where  Vi  Is  the  value  of  the  property  P at 
the  temperature  Ti  (1  • 1,  2,  3,  ••:) 

Option  2 

HI  » number  of  polynomial  coefficients  di  to  be  Input  by  the  user  as 
an  alternate  to  Inputting  (Ti,  Vi)  pairs 

ai  ■ coefficients  of  a polynomial  equation  which  gives  the  value 
of  the  property  P as  a function  of  temperature  (i  ■ 0,  1, 

2,  3,  •••) 

Currently,  there  are  ten  properties  which  must  be  Input  for  each  material 
using  one  of  the  two  options  described.  Assigned  property  numbers  P 
and  the  corresponding  property  descriptions  are  given  In  Table  5. 


TABLE  5.  MATERIAL  PROPERTY  IDENTIFICATION 


p 

Property  Description 

Symbol 

1 

Elastic  modulus  (Young's) 

E 

2 

Secant  modulus  at  point  A 

Ea 

3 

True  stress  at  point  A 

dA 

4 

True  stress  at  rupture 

*r 

5 

True  strain  at  rupture 

€r 

6 

Poisson's  ratio 

V 

7 

Coefficient  of  thermal  expansion 

*T 

8 

Mass  density  (wt  per  unit  vol.  x 1/acc.  of  grav.) 

P 

9 

Specific  heat 

s 

10 

Damping  to  stiffness  ratio 

h 

The  first  five  properties  define  the  Ramberg-Osgood  representation 
of  the  true  stress-true  strain  curve,  reference  Figure  9.  Current 
formulations  are  based  on  the  assumption  that  the  materials  are 


Isotropic.  Sources  for  material  properties  are  References  3 and  4. 
See  also  pages  29  through  31. 

' 
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Data  Code  11,  Bird  Impact  Definition  Data 

User  inputs  cards  as  described  below  and  one  end  card.  See  pages  81 
through  88  for  discussion,  definitions  and  equations. 


RT  ADJ  INTEGER 

FLOATING  POINT 

a 

a 

a 

a 

c 

d 

e 

12 

3 

7 

n 

15 

19 

34 

49 

E 

9999 

| Er 

d Ca 

rd 

Card  1 


a - 1 

b * row  dimension  for  output  matrix  of  incremental  loads 
3 3 x number  of  joints  in  the  model 
c * bird  mass  m 
d * bird  velocity  v 

e 3 effective  length  L » reference  Figure  12 
Card  2 
a 3 2 

b 3 number  of  time  increments  n^  during  which  impact  is  applied 
c,  d,  e 3 X,  Y,  Z components  of  unit  vector  b parallel  to  the  bird  path 
Card  3 
a 3 3 

b 3 not  used 

c,  d,  e ■ X,  Y,  Z components  of  a unit  vector  IT  normal  to  the  surface 
of  the  structure  at  the  impact  point 
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Card  4 


a ■ 4 

b ■ not  used 

c,  d,  e « X,  Y,  Z components  of  a unit  vector  T in  the  same  plane  as  "b 
and  n^  in  the  approximate  direction  of  motion  of  the  Impact 
footprint,  and  normal  to  the  desired  direction  of  the  re- 
sultant impact  load,  reference  Figure  12. 

Card  5 
a * 5 

b ■ not  used  • 

c,  d,  e * X,  Y,  Z coordinates  of  the  impact  point 
Cards  6 (There  are  n^  of  these  cards) 
a » 6 

b a time  increment  identification  number  0 (/3  * 1,  2,  3 •••,  Njg) 

c » distance  Dg  traveled  by  the  Impact  footprint  in  the  time 
period  from  tQ  to  tg_j 

d ■ load/time  distribution  factor  fg  for  Increment  8 
e ■ not  used 
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Data  Code  12,  Impact  Footprint  Joint  List 


User  inputs  nj  of  these  cards  and  one  end  card,  see  page  87. 


RIGHT  ADJUSTED  INTEGER  DATA 


End  Card 


Card  col  72  non-blank  for  continuation 


j3  » time  increment  identification  number  ( j3  ■ 1,  2,  3,  •••,  n j ^ ) 


j-jp  * joint  numbers  of  all  joints  to  which  the  resultant  bird  impact 
force  is  to  be  distributed  in  the  pth  time  increment 


Data  Code  20,  Axial  Bar  Element  Definition 


User  Inputs  one  card  for  each  bar  and  one  end  card.  See  pages  5 
through  10  for  discussion  of  finite  elements. 


N ■ bar  Identification  number 

p,  q ■ joint  numbers  definlna  the  ends  of  the  bar 

M ■ material  Identification  number  (references  code  10  data) 

A ■ cross-sectional  area  of  the  bar 


A 


Data  Code  30,  Membrane  Element  Definition 

User  Inputs  one  card  for  each  membrane  and  one  end  card.  See  pages 
5 through  10  for  discussion  of  finite  elements. 


N * membrane  Identification  number 

p,  q,  r,  s * joint  numbers  defining  the  four  comers  of  the 
membrane,  see  below 


M * material  Identification  number  (references  code  10  data) 

t » thickness  of  membrane 

5pq  * stress  orientation  angle  in  degrees 


Data  Code  40,  Cell  Element  Definition 

User  Inputs  one  card  for  each  cell  and  one  end  card  . See  pages 
5 through  10  for  discussion  of  finite  elements. 


Data  Code  50,  Point  Mass  Element  Definition 


User  inputs  one  card  for  each  point  mass  and  one  end  card.  See  page 
10  for  discussion. 


N * point  mass  identification  number 

p = joint  number  defining  the  location  of  the  point  mass 

Mpx»  MpY , MpZ  = uncoupled  mass  contributions  which  will  be  added 
to  the  diagonal  of  the  structural  mass  matrix  in  the  X,  Y,  Z 
degrees  of  freedom  at  joint  p 
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FORMAT  Phase  2 Data 


Three  data  deck  setups  are  required  as  indicated  in  Figure  10. 

It  is  not  necessary  for  the  user  to  refer  to  the  FORMAT  documentation 
to  use  these  setups.  In  keypunching  the  cards,  it  is  important  to  begin 
in  the  card  columns  shown.  Cards  beginning  with  the  % character  in  CC1 
have  a control  card  name  field  In  CC2-15  and  an  option  field  in  CC16-72. 
The  control  card  name  and  the  optional  information  must  be  entered 
left  justified  In  these  fields  as  shown.  All  other  cards  in  the 
sequence  must  begin  In  card  column  7 and  may  be  punched  in  free  form. 

An  exception  is  matrix  data  which  has  the  following  card  format. 


Header  Card  (one  required  for  each  matrix) 


G 

IS 

Name 

8 

up 

m 

Data  Card  (use  as  many  as  required  to  input  non-zero  elements  of  matrix) 


12  - 5b  - 9 


Value 


Exp 


h 1-1 920-221 


24-2728-311 


Value  Exp 


33-41W2-44 


46-49150-53 


Value  Exp  Name 


55-63B4-66b 


End  of  Matrix  Data  Card  (one  required) 


Notes:  1.  "Value"  Is  of  the  form  tXXX.XX,  plus  sign  may  be  omitted. 

2.  "Exp"  Is  of  the  form  ±XX,  "Exp"  * ±00  may  be  omitted. 

3.  More  than  one  matrix  may  be  Input,  but  only  one  "E"  card 
is  required  after  the  last  matrix. 
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The  following  definitions  apply. 


i 


max,  '■'max 


row  and  column  dimensions  of  the  matrix 


Name  = one  to  six  character  alphameric  matrix  name  assigned 

by  the  user.  The  first  character  must  be  alphabetical. 
Enter  free  form. 


1,  j * right  adjusted  integer  values  specifying  the  row  and 

column  location  of  a matrix  element  whose  value  is 
given  in  the  next  field 


Value,  Exp 


value  of  the  i,  j matrix  element 


FORMAT  Deck  Setup  1 


Card  Card  Column 
No.  1 7 16 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


SF0RMAT  STANDARD 

$RUN  G0.L0GIC 

INPUT  TAPE  (MTAPE.l) 

INPUT  TAPE  (LTAPE.l) 

OUTPUT  TAPE  (DEC0MP.1) 

S I ZE ( NW0RK , K0NST ) 

SINSTRUCTION 

PT.DPT  = DP0T  -DEJC0L-  (1) 

SAVE(DEC0MP)DP0T 

MR  = PRUF.MEL  -SEQWF. 

PRUPJ, PRUPS  = PRUPT  -DEJCOL-  (3r>j) 
PR  = PRUPJ  -MULT-  PT 

PRT  = PR  -TRANSP. 

X.LTL.UR  = PRUF.KEL  -SEQWF-  PRT 
UT  = PRUP  -TMULT-  UR 

SAVE( DEC0MP) MR, PRUPJ ,LTL  ,UR 
PRINT ( , , ,)UT 

SEND 


An  alternate  setup  after  card  17  which  will  accommodate  user  input 
matrices  on  cards  Is: 

SMATRIX  LIST 

Card  Input  matrix  data,  see  pages  65  and  66  for  card  formats. 

SEND 

In  card  6 NW0RK  represents  words  of  working  core  storage,  and  K0NST 
Is  the  maximum  matrix  dimension  allowed.  Default  values  are  10000  and 
2000.  Estimates  can  be  calculated  from  equations  12  and  13,  where  w 
Is  the  wave  front,  see  "Joint  Numbering"  on  page  28,  and  npp  Is  given 
by  equation  22. 

In  card  11 , 3nj  ■ 3 x number  of  joints  In  the  model. 
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NWORK  * w2/2 
KONST  i nQC 


FORMAT  Deck  Setup  2 


Card  Card  Colurm 
No.  1 7 


$F0RMAT 

$RUN 


STANDARD 


SRUN  G0,L0GIC 

INPUT  TAPE  (MTAPE.1) 

INPUT  TAPE  (DEC0MP  ,1) 

0UTPUT  TAPE  (TRANSF.l ) 

S IZE ( NW0  RK , K0NST ) 

SINSTRUCTI0N 

VAl.TR  = MR  -USER06-  LTL 
PRINT ( ,,  ,)VAL 
PBUPJ  = TR  -TMULT-  PRUPJ 
PBUF  = TR  -TMULT-  PRUF 
DPB0  = PBUPJ  -MULT-  DP0T 
SAVE(TRANSF)TR , PBUPJ, PBUF, DPB0 
M0DES  > PRUPJ  -TMULT-  TR 
PRINT( ,,F2 ,.O5)M0DES 
SSPECIAL 


CC6  CC12 


CC24 

1.0E-6 


(Rt.  adj.  to  noted  CC) 


In  card  15,  the  value  .05  Is  a print  cutoff  which  provides  that 
any  matrix  element  whose  value  Is  - .05  will  not  be  printed.  The 
user  may  change  the  cutoff  value  If  desired.  If  the  value  Is  zero  or 
blank,  all  matrix  elements  will  be  printed  except  those  whose  value 
Is  zero. 


1 


In  card  17,  nQ  is  the  number  of  modes  to  be  used  in  the  modal 
transformation.  The  modes  extracted  will  correspond  to  the  n^  lowest 
natural  frequencies  of  the  structure.  The  value  1.0E-6  is  the  con- 
vergence criterion  for  the  eigenvalue  calculation  and  is  based  on 
current  experience  levels-. 


FORMAT  Deck  Setup  3 


Card  Card  Coluim 
No.  1 7 


16 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


$ FORMAT  STANDARD 

$RUN  G0 .LOGIC 

INPUT  TAPE  (MTAPE.l) 

INPUT  TAPE  (TRANSF.l) 

OUTPUT  TAPE  (FILE20 ,1 ) 

OUTPUT  TAPE  ( FILE21 ,1) 
SIZE(NW0RK,K0NST) 

SINSTRUCTION 

SAVE(FILE20)PBUF 
SAVE(FIL£20)MPT 
SAVE ( FILE20 ) UZER0 

SAVE( FILE20 )ECT ,MEL ,KBEL ,CBEL ,FFBAR 
SA VE ( FILE20 ) S I GFB .EPS I G .DEBT , E VT , C0NST 
SAVE( FILE21  )DPB0 
SAVE ( FILE21 )PBUPJ 

SEND 


Incremental  Response  Data 

Card  input  to  the  Incremental  Solution  program  consists  of  one 
card  to  specify  options  and  dimensional  quantities  followed  by  one 
or  more  cards  as  needed  to  specify  the  time  points  which  mark  the 
end  of  each  increment. 
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Card  1 (card  colutms  1 through  37  are  used) 


r 


RIGHT  ADJUSTED  INTEGER  DATA 


s 

nEL 

% 

nG 

I 

n0 

UW0RK 

1 - 6 

7 - 12 

13  - 18 

19-24  j 

1 

26  - 31 

32  - 37 

S 3 starting  increment  number 

3 1 

n^L  3 total  number  of  elements  in  the  model  including  bars,  mem- 
branes, cells,  point  mass. 
n3  3 nunfcer  of  increments  in  the  solution 
n^  3 number  of  transformation  modes 

T,  F 3 true,  false,  T means  that  element  temperature  changes  will  be 
calculated  as  a result  of  damping  energy  converted  to  heat. 

F means  that  they  will  not  be  calculated. 

3 total  number  of  joints  in  the  model. 

JW0RK  3 size  of  core  work  space  required  and  must  satisfy  the  following: 

JW0RK  > 1200  + 10.5  n*  + 29.5  nQ  + "EL  + ne  + 3 nj 
JW0RK  £ 7000  + 2.5  n^  + 115.5  nG  + nEL  + ns  + 4 nJ 
Card(s)  2 (card  columns  1 through  72  are  used) 


FLOATING  POINT  DATA 


tl 

h 

• 

• 

t8 

1 - 12 

13  - 24 

25  - 36 

37  - 48 

49  - 60 

61  - 72 

t.  » time  from  beginning  of  impact  (trt)  to  the  end  of  the  eth  incre- 

p Q 

ment,  s 3 l,2,3,...,na  . The  t.  are  the  time  points  at  which 

P P 

response  solutions  are  calculated. 


The  t must  be  input  in  sequence  at  the  rate  of  six  values  per  card 

P 

until  ng  values  have  been  entered. 

' 
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Postprocessor  Data 

Card  input  to  the  Postprocessor  program  consists  of  six  option 
definition  cards,  four  title  cards,  and  a variable  number  of  selectivity 
cards. 


Card  Format  for  Cards  1 through  6 


RIGHT  ADJUSTED  INTEGER  DATA 

FLOATING  POINT 

RIGHT  ADJUSTED  ; 

a 

x 

f 

c 

1 

e 

S3 

9 

h 

D 

a 

3 - 10 

B 

E 

E 

E 

B 

16  - 20 

21-30 

31  - 40 

-i 

INTEGER  Dm. A 

c — 

i 

j 

k 

X 

i— 

41  - 50 

51  - 60 

61  - 70 

71  - 80 

Card  1.  Selectivity  Options 
a 3 1 

b 3 0 or  1.  "Zero"  means  that  the  user  desires  printout  for 

all  time  increments.  "One"  means  that  the  user  will  specify 
selected  time  increments  for  which  he  desires  printout, 
c 3 0 or  1 . "Zero"  means  that  the  user  desires  printout  for  all 
joints  as  dictated  by  card  2.  "One"  means  that  the  user  will 
specify  selected  joints  for  which  he  desires  printout, 
d 3 0 or  1 . "Zero"means  that  the  user  desires  printout  for  all 
bars  as  dictated  by  card  4.  "One"  means  he  will  specify 
selected  bars. 

e 3 0 or  1.  "Zero"  means  that  the  user  desires  printout  for  all 
membranes  as  dictated  by  card  5.  "One"  means  he  will  specify 
selected  membranes. 
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f 


I 


g 


h 


i 


j 


k 


0 or  1 . "Zero"  means  that  the  user  desires  printout  for  all 
cells  as  dictated  by  card  6.  "One"  means  he  will  specify 
selected  cells. 

value  which  the  program  will  use  to  scale  the  joint  displace- 
ments. May  be  useful  when  plotting  the  deformed  structure 
joint  coordinates.  If  this  field  is  blank,  the  default  value 

for  the  scalar  is  1.0.  The  value  is  determined  by  experience. 

KW0RK,  which  is  the  size  of  work  space.  The  default  is  i 0000. 
10000  < KW0RK  > 640  + L]  + L2  + 2ng  + 3n^  + 3nfi 

where  6ng  < >60  and  l*2  is  the  larger  of  the  numbers  of 

joints,  bars,  membranes  or  cells. 

KLNMAX,  which  is  the  maximum  number  of  lines  per  page, 
measured  from  the  top  of  the  page.  The  minimum  value  is  20, 
and  a blank  field  results  in  a default  value  of  60. 

KCLMAX , which  is  the  number  of  elements  across  the  page.  The 
minimum  value  is  1,  and  a blank  field  results  in  a default 
(maximum)  value  of  8. 

IFLGPR,  a non-blank  field  sets  a checkout  print  flag  for 
intermediate  print- -normally  not  turned  on  by  the  user. 


Card  2.  Joint  Data  Options 
a = 2 

b,c,d,e  are  individually  set  equal  to  0 or  1 to  control  printing 

of  joint  coordinates,  displacements,  velocities,  and  acceler- 
ations, respectively.  "Zero"  means  no  print.  "One"  means 
print. 

f through  k are  not  used. 

Card  3.  Modal  Response  Options 
a a 3 

b,c,d  are  individaully  set  equal  to  0 or  1 to  control  printing  of 

modal  displacements,  velocities  and  accelerations,  respectively. 
"Zero"  means  no  print.  "One"  means  print. 

< through  k are  not  used. 
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Cards  4,  5,  6.  Element  Data  Options 


a = 4 for  bar  elements. 

* 5 for  membrane  elements. 

= 6 for  cell  elements 

b,c,d,e  are  individually  set  equal  to  0 or  1 to  control  printing 

of  element  forces,  stresses,  strains,  and  equivalent  stresses, 
respectively.  "Zero"  means  no  print.  "One"  means  print, 
f through  k are  not  used. 

Card  format  for  cards  7,  8,  9 and  10. 


RT  ADJ  INTEGER 

ALPHAMERIC 

a 

b 

D 

a 

3 - 10 

11  - 72 

L 

These  four  cards  allow  the  user  to  specify  two  120  character  lines 
of  title  information  to  be  printed  at  the  top  of  each  page  of  printout. 

a = 7,  b = characters  1-62  of  first  title  line, 

a = 8,  b = characters  63-124  of  first  title  line, 

a * 9,  b = characters  1-62  of  second  title  line, 

a = 10,  b - characters  63-124  of  second  title  line. 

Card  format  for  selectivity  cards. 


Increment  selections--requi red  only  when  card  1 has  a value  of  1 
in  column  11 . 


a 


= 20 

= list  of  time  increment  numbers  3 for  which  printout  is  desired. 
Use  additional  cards  as  necessary, 
b = 9999  in  the  last  card. 

Joint  selections — required  only  when  card  1 has  a value  of  1 in 
column  12. 

a = 30 

c^  * list  of  joint  numbers  for  which  printout  is  desired.  Use 
additional  cards  as  necessary, 
b = 9999  in  the  last  card. 

Bar  selections  — requi red  only  when  card  1 has  a value  of  1 in 
column  13. 

a = 40 

c^  = list  of  bar  numbers  for  which  printout  is  desired.  Use 
addi ti onal' cards  as  necessary, 
b = 9999  in  the  last  card. 

Membrane  se1ections--required  only  when  card  1 has  a value  of  1 
in  column  14. 

a a 50 

c.|  = list  of  membrane  numbers  for  which  printout  is  desired.  Use 

additional  cards  as  necessary, 
b = 9999  in  the  last  card. 

Cell  selections — required  only  when  card  1 has  a value  of  1 in 
colurm  15. 

a * 60 

c^  * list  of  cell  numbers  for  which  printout  is  desired.  Use 
additional  cards  as  necessary, 
b * 9999  in  the  last  card. 
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In  any  of  the  selectivity  cards  described  above,  an  inclusive  range 
of  numbers  may  be  specified  by  entering  the  lower  and  upper  bounding 
numbers  with  a non -blank  character  (such  as  — ) entered  between  them  in 
the  first  card  column  of  the  field  containing  the  upper  bounding  number. 

Job  Control  Language  (JCL) 

The  Programing  Manual,  which  is  Part  3 (a  separate  volume)  of  this 
report,  describes  program  implemenation  procedure  at  a given  computer 
installation  and  the  job  control  language  (JCL)  required  to  execute 
the  IMPACT  computer  programs.  The  JCL  deck  contains  estimates  of  computer 
time,  core,  and  other  resource  requirements  for  the  job,  manipulates 
input  and  output  data  files,  and  calls  and  executes  program  load  modules. 

A JCL  deck  together  with  the  appropriate  input  data  cards  must  be  set 
up  and  loaded  into  the  coirputer  to  perform  the  functions  shown  in 
Figure  10.  Depending  on  the  size  of  the  finite  element  model  being 
analyzed,  the  JCL  setup  may  be  tailored  to  perform  each  function  as  an 
individual  computer  run  in  the  case  of  a large  model  or  to  perform  all 
functions  in  one  computer  run  for  small  models.  Reference  Appendix  A 
for  additional  Information  and  examples  of  JCL. 

INPUT  PREPROCESSING 

Input  preprocessing  includes  the  input  and  editing  of  model  defini- 
tion decks  by  means  of  the  CDC  UPDATE  procedure,  the  generation  of 
additional  model  definition  data  using  the  Laminate  Generator,  and  the 
generation  of  the  bird  impact  incremental  loads  matrix  using  the  Loads 
Generator  program.  The  makeup  of  the  model  definition  deck  is  dis- 
cussed under  Data  Code  Descriptions  in  the  previous  subsection.  UPDATE 
is  a CDC  system  module  for  creating,  accessing  and  editing  data  files 
and  Is  utilized  through  appropriate  JCL. 

Laminate  Generator 

The  use  of  the  Laminate  Generator  program  Is  not  manditory,  but  it 
is  Intended  to  relieve  the  user  of  much  of  the  effort  required  to  pre- 
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pare  model  definition  data  for  the  transparency,  especially  when  it  is 
a multi-layered  laminate  having  compound  curvature. 


General  Description 

In  order  to  take  advantage  of  the  Laminate  Generator  capability 
the  transparency  must  be  modelled  entirely  of  cell  elements.  Minimum 
required  input  data  includes  data  code  2 joint  coordinates  of  the  joints 
in  the  surface  of  the  transparency,  data  code  7 laminate  definition 
cards  and  data  code  40  cards  giving  the  cell  identification  number 
and  the  first  four  corner  joint  numbers  for  each  cell  in  the  surface 
layer  of  the  transparency.  The  program  then  calculates  joint  coor- 
dinates for  all  additional  joints  required  to  define  the  transparency 
geometry  and  produces  cell  definition  data  to  complete  surface  cell 
definitions  and  to  define  all  additional  cell  elements  needed  to  form 
sublayers  of  the  transparency. 

Consider  a 5-layered  laminated  transparency  modelled  with  a 10  x 10 
surface  grid.  Such  a model  would  have  600  joints  and  405  cell  elements. 
The  user  would  input  joint  coordinates  for  100  surface  joints,  partial 
cell  definitons  for  81  cells  in  the  surface  layer,  and  information 
specifying  thicknesses  and  materials  of  layers  and  numbering  of  joints 
and  cells.  The  program  will  produce  the  completed  data  code  2 deck 
for  600  joints,  a data  code  6 deck  defining  unit  vectors  normal  to  the 
surface  at  the  100  surface  joints,  and  a completed  data  code  40  deck 
defining  the  405  cell  elements. 

The  input  file  is  Fortran  logical  unit  1 and  is  passed  to  the 
Laminate  Generator  from  an  UPDATE  data  file,  as  indicated  in  Figure  10. 
Table  4 shows  optional  data  code  decks  which  may  be  present  in  the 
input  file  as  well  as  the  minimum  required  input  discussed  above.  If 
the  optional  input  is  present,  it  will  be  duplicated  (except  data  codes 
7 and  3)  in  the  output  file,  thus  creating  a complete  model  definition 
deck.  This  completed  model  definition  deck  is  now  in  the  form  required 
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for  use  in  the  Loads  or  Initial  Generators.  If  the  optional  input  is 
not  present,  the  Laminate  Generator  output  data  file  will  have  to  be 
merged  with  other  data  files  as  appropriate  to  form  a complete  model 
definition  file.  The  output  file  is  Fortran  logical  unit  2 and  is 
copied  to  an  UPDATE  data  file  and  saved  for  future  use  or  it  may  be 
copied  to  appropriate  Fortran  logical  units  and  passed  directly  to  the 
Loads  or  Initial  Generators. 

The  Laminate  Generator  unconditionally  prints  out  a listing  of  all 
data  code  decks  present  in  the  output  file. 


Definitions  and  Use  Details 

Figure  11  shows  a representative  portion  of  a finite  element  model 
for  a 3-layered  laminate  and  support  structure.  The  following  defin- 
itions apply. 

1)  Surface  normal --a  unit  vector  directed  into  the  laminate  based  on 
the  average  of  cross  products  of  surface  element  edge  vectors 
emanating  from  a surface  joint. 

2)  Joint  stack — a set  of  joints  including  a surface  joint  and  all 
joints  through  the  laminate  to  the  opposite  surface,  reference 
Figure  11  (1 ,2,3,4),  (13,14,15,16),  etc.  The  joint  stack  line  is 
usually  parallel  to  the  surface  normal  as  defined  above. 

3)  Layer  thicknesses--thi ckness  t input  by  the  user  for  each  layer  i 
and  used  by  the  program  together  with  normal  vectors  and  surface 
joint  coordinates  to  calculate  the  coordinates  of  all  other  joints 
In  each  stack. 

4)  Joint  numbering  1ncrement--an  Integer  value  aj^  input  by  the  user 
for  each  layer  i and  used  by  the  program  to  assign  joint  numbers 
in  the  joint  stacks  starting  from  the  user  assigned  surface  joint 
number  and  proceeding  through  the  stack.  In  Figure  11,  aJ^  3 1 for 
all  layers.  The  value  for  aJz  cannot  be  negative. 
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5)  Element  stack— a set  of  cell  elements  including  a surface  layer 
cell  and  a cell  for  each  additional  layer  in  the  laminate,  reference 
Figure  11,  (1,7,13),  (6,12,18),  etc. 

6)  Element  numbering  increment— an  integer  value  aE^  Input  by  the  user 
for  each  layer  £ and  used  by  the  program  to  assign  element  numbers 
in  the  joint  stacks  starting  from  the  user  assigned  surface  element 
number  and  proceeding  through  the  stack.  In  Figure  11,  aE^  * 6 for 
all  layers.  The  value  for  ^ cannot  be  negative. 

7)  Surface  offset— a distance  tQ  input  by  the  user  and  used  by  the 
program  to  calculate  new  surface  joint  coordinates  such  that  the 
entire  surface  is  moved  a constant  distance  tQ  outward  or  Inward 
from  its  original  position  (positive  t0  results  in  offset  opposite 
to  the  unit  normal  vector  directions). 

8)  Partial  cell  element  definitions— a code  40  card  input  by  the  user 
having  the  cell  identification  number,  the  first  four  comer  joint 
numbers,  and  the  stress  orientation  angle  (if  applicable)  entered 
In  their  appropriate  fields.  Note  that  the  comer  joint  numbers 
must  be  surface  joints  and  they  must  be  entered  in  a clockwise 
sequence  looking  down  on  the  surface  in  the  direction  of  the 
unit  vectors.  In  Figure  11  for  cell  number’ 1 the  sequence  could  be 
one  of  four  choices:  (1  ,24,28,5),  (24,28,5,1),  (28,5,1  ,24)  or 
(5,1  ,24,28). 

9)  Variable  layer  thlcknesses-thickness  tt  input  by  the  user  for 
each  layer  i at  each  surface  joint  as  necessary  to  define  a 
realistic  variation  In  layer  thickness. 

Items  defined  in  3,  4,  6 and  7 above  are  Input  on  data  code  7 
cards  together  with  a material  reference  for  each  layer.  The  variable 
thickness  data  defined  in  9 is  input  on  data  code  8 cards. 

As  mentioned  in  2 above,  the  joint  stack  line  will  be  parallel  to 


i 
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the  surface  normal  vector  if  only  the  laminate  surface  joint  coordinates 
are  input.  However,  the  user  may  control  the  direction  of  the  stack 
line  by  inputting  coordinates  for  the  last  joint  in  the  stack  as  well 
as  the  surface  joint.  In  fact,  intermediate  joint  coordinates  may  also 
be  input  if  desired.  In  Figure  11  this  Is  illustrated  along  the  edge  of 
the  laminate  which  is  attached  to  the  jpport  structure  at  joints  16, 

39  and  62.  User  input  of  coordinates  for  these  joints  controls  the 
stack  line  direction  to  insure  that  the  attach  line  between  the  lami- 
nate and  the  support  structure  will  be  where  the  user  wants  it.  The 
undesirable  alternative  would  be  to  adjust  coordinates  of  surface 
joints  13,  36  and  59  to  force  the  surface  normals  to  intersect  the 
desired  attach  line. 

The  Laminate  Generator  discussion  is  summarized  by  again  referring 
to  Figure  11.  For  this  model  the  user  should  input  data  code  2 coordi- 
nates for  laminate  surface  joints  1 , 5 , 9,  13  , 24  , 28  , 32  , 36  , 47  , 51  , 

55  and  59  and  for  support  structure  joints  16-23,  39-46,  and  62-69. 

He  should  input  data  code  7 to  define  surface  offset  (if  necessary), 
layer  material,  joint  and  element  numbering  increments,  and  thickness 
for  each  of  the  three  layers.  Data  code  8 cards  should  be  input  if 
there  are  variations  in  thickness  of  the  layers.  Finally,  data  code 

40  cards  should  be  input  to  partially  define  surface  elements  1-6  and 
completely  define  elements  19-24.  Optional  data  codes,  1,  3,  4,  5, 

10,  11  , 12,  20,  30,  40,  and  50  should  be  included  if  available.  Now, 
execution  of  the  Laminate  Generator  program  will  result  in  a completed 
data  code  2 deck  with  missing  joint  coordinates  inserted,  a new  data 
code  6 deck  of  surface  normal  vectors  and  associated  surface  areas, 
and  a completed  data  code  40  deck  with  all  missing  element  definition 
data  Inserted.  In  addition,  all  optional  data  code  decks  which  were 
present  in  the  input  file  will  have  been  output  in  proper  sequence. 

Data  codes  7 and  8 which  are  used  by  the  Laminate  Generator  are  excluded 
from  the  output. 
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Loads  Generator 


The  problem  of  determining  a space/time  distribution  of  applied 
loads  at  the  model  joints  which  accurately  represents  bird  impact 
phenomena  is  a very  difficult  one.  Reference  2 reports  on  experimental 
and  analytical  investigations  of  bird  impact,  presents  some  practical 
formulations  of  force/time  distribution  as  well  as  pressure  distri- 
butions within  the  impact  footprint,  and  makes  recommendations  for 
the  continuation  of  this  important  work. 


General  Description 

The  Loads  Generator  program  is  an  elementary  capability  which 
calculates  the  average  impact  force  using  the  formulation  of  Reference 
2,  allows  the  user  to  specify  the  force/time  distribution  and  the 
movement  of  the  footprint,  and  calculates  the  resulting  load  com- 
ponents on  joints  within  the  footprint.  The  load  distribution  to  the 
joints  is  not  based  on  impact  pressures  but  instead  is  based  on  the 
method  of  Lagrangian  multipliers  to  solve  equilibrium  equations  relat- 
ing joint  loads  to  the  resultant  impact  force.  Target  compliance  is 
not  accounted  for,  but  a travelling  footprint  can  be  acconmodated. 

The  input  file  is  Fortran  logical  unit  1 and  is  passed  to  the 
Loads  Generator  from  an  UPOATE  data  file,  as  indicated  in  Figure  10. 
Data  code  decks  required  to  be  present  in  the  input  file  are  shown 
in  Table  4.  Optional  data  codes  present  in  the  input  will  be  read 
over  and  Ignored.  The  output  file  is  Fortran  logical  unit  3 and 
contains  the  incremental  loads  matrix  stored  in  the  form  of  a FORMAT 
matrix  tape  named  LTAPE.l.  It  may  be  copied  either  to  the  appropriate 
Fortran  logical  unit  for  input  to  a F0RMAT  step  or  to  tape  and  saved 
for  future  use.  Alternately,  the  output  file  which  is  also  stored  in 
card  image  form  on  Fortran  logical  unit  2 may  be  punched  on  cards  in 
a form  suitable  for  matrix  card  input  to  the  FORMAT  step.  These 
options  are  exercised  by  means  of  JCL  cards. 
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The  Loads  Generator  unconditionally  prints  a load  summary  for 
each  time  increment  during  the  impact.  The  suimary  includes  the 
load  (time)  increment  number,  bird  mass  and  velocity,  duration  of 
impact  and  average  impact  force,  direction  and  magnitude  of  the 
resultant  impact  force  for  the  given  increment,  the  loaded  joint 
numbers  together  with  their  total  loads  for  the  given  increment  and 
their  incremental  load  changes  from  the  previous  increment,  and 
finally  a column  of  the  output  Incremental  loads  matrix. 

Definitions  and  Use  Details 

Figure  12  shows  parameters  and  diagrams  which  the  user  must 
consider  in  order  to  prepare  input  data  for  the  Loads  Generator.  The 
following  definitions  apply. 

1)  3nj  — row  dimension  of  incremental  loads  matrix 

2)  m — bird  mass 

3)  v — bird  velocity 

4)  L — estimated  total  travel  of  footprint  along  d from  point  of 

impact  to  point  of  detachment  divided  by  cosine  9 

5)  9 — angle  through  which  the  bird  material  is  expected  to  be 

deflected  during  the  impact 

6)  b — unit  vector  parallel  to  bird  path  prior  to  impact 

7)  n — unit  vector  normal  to  transparency  surface  at  Impact  point 

8)  d — unit  vector  In  plane  with  vectors  b and  n and  making  the 

angle  e with  vector  i) 
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9)  S — time  increment  identification  number 

10)  tg  — time  at.  the  end  of  the  sth  time  increment 

11)  Fg  — resultant  impact  force  at  the  footprint  c.g.  acting  during 

the  time  increment 

12)  F — average  impact  force 

13)  fg  --  non-dimensional  impact  force  factor  for  the  increment 

equal  to  the  ratio  Fa/F 

p avy 

14)  Dg  --  distance  from  the  impact  point  to  the  footprint  c.g.  at 

the  beginning  of  the  6th  increment  measured  parallel  to 

vector  d (for  8 * 1 , D-j  * 0) 

15)  n^3  — number  of  time  increments  in  the  impact  duration  time 

16)  t^  --  impact  duration  time 

17)  XQ,  Y , ZQ  --  coordinates  of  the  impact  point 

Data  codes  11  and  12  are  used  to  identify  bird  impact  input  data 
for  the  Loads  Generator.  The  following  discussion  is  intended  to 
clarify  the  details  of  preparing  this  input. 

The  value  njj  is  the  total  number  of  joint  translational  degrees 
of  freedom  in  the  finite  element  model  and  is  found  from: 


nTJ  3 3nJ 

where  nj  is  the  total  number  of  joints  in  the  model. 


04) 


Bird  mass  m and  velocity  v are  specified  by  the  user  in  order 
to  produce  the  desired  impact.  It  should  be  recognized  that  v is 
the  relative  velocity  of  bird  to  transparency. 
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The  value  of  L is  an  important  parameter  because  it  determines 
the  impact  duration  time.  It  should  be  estimated  based  on  the  user's 
experience  and  best  judgement  for  each  impact  case.  Reference  2 defines 
L as  the  length  of  the  bird  in  the  case  of  a rigid  target  whose  surface 
is  normal  to  the  bird  path,  and  gives  a modified  effective  length  in 
the  case  of  inclined  surfaces.  The  travelling  footprint  concept 
described  here  is  a further  modification. 


At  this  point  the  user  must  establish  the  number  of  impact  incre- 
ments nj3  and  the  time  points  tg  marking  the  end  of  each  increment. 
Figure  13  is  a nondimensional  representation  of  the  triangular  impact 
force/time  distribution  as  derived  in  Reference  2.  It  shows  the  impact 
divided  into  ten  equal  increments  and  gives  the  values  of  the  force 
factor  fa.  It  is  not  necessary  that  the  increments  be  equal.  The 

P 

values  for  nT_,  ta  and  f.  should  be  determined  in  as  much  detail  as 

• Pm  P 

necessary  to  reasonably  represent  the  desired  force/time  distribution, 

whether  it  be  triangular  or  some  other  shape.  For  example,  the  dotted 

line  in  Figure  13  represents  a trapezoidal  distribution  that  may  be 

more  applicable  for  the  compliant  target  and  travelling  footprint. 

The  value  of  tA  is  needed  to  determine  the  actual  values  of  ta. 

d S 


V7 


05) 


To  insure  that  the  correct  total  impulse  is  being  input  to  the 
Loads  Generator  the  user  should  check  that  the  following  equation  is 
satisfied  (s 


1 ,2,3,...  ,njg). 


ZVWi> 


(16) 


The  unit  vector  b,  n and  d are  input  in  terms  of  their  x,  y,  z 
components  in  the  model  coordinate  axes  system.  The  user  must  calcu- 
late the  components  of  b consistent  with  the  desired  bird  path.  The 
components  of  n can  be  determined  from  the  data  code  6 surface  normals 
for  the  joints  nearest  to  the  point  of  impact.  The  data  code  6 infor- 


85 


w'. 


AD-A062  61t 


UNCLASSIFIED 


2*2 
ADA 
0626 1 A 


DOUGLAS  AIRCRAFT  CO  LONG  BEACH  CALIF  F/G  5/2 

AIRCRAFT  WINDSHIELD  BIRO  IMPACT  MATH  MODEL.  PART  2.  USER'S  MANU— ETC(U) 
DEC  77  G R EIDE  F33615-75-C-3105 


MDC-J7174-PT-2 


AFFDL-TR-77-99-PT-2 


END 


DATE 

FILMED 


nation  is  part  of  the  unconditional  printout  from  the  Laminate  Generator. 
The  components  of  d are  given  by  the  following  equation. 

d « (cos  8)  b - (sine)  ^ j *-|j  (17) 

The  value  for  8 can  be  scaled  from  a layout  of  the  bird  impact  and  trans- 
parency geometry  similar  to  Figure  12. 

The  coordinates  XQ,  Y0,ZQof  the  impact  point  are  in  the  model 
coordinate  axes  system  and  represent  the  intersection  of  the  bird  path 
line  and  the  transparency  surface. 

The  following  equation  gives  the  values  of  D.  for  8 = 1 ,2,3,.. ,fnr 

D IS. 

°3  = Vl  v C0S  0 08) 

In  order  to  determine  the  joint  numbers  which  lie  within  the  foot- 
print boundary  for  each  increment  the  user  should  lay  out  a view  of 
the  transparency  surface  showing  the  impact  point,  the  model  grid  and 
surface  joint  numbers  for  that  portion  of  the  surface  over  which  the 
bird  material  will  pass,  see  Figure  12.  The  intersection  of  the  plane 
containing  vectors  b,  n and  d with  the  surface  should  be  drawn  on 
the  layout,  and  the  distances  0o  should  be  marked  off  from  the  impact 
point  along  the  intersection  line  to  locate  the  footprint  c.g.  for 
each  increment.  The  footprint  boundary  ellipse  should  be  drawn  around 
each  c.g.  point  Including  the  impact  point.  The  minor  axis  of  the 
ellipse  may  be  assumed  to  be  equal  to  the  diameter  of  the  bird.  The 
major  axis  Is  equal  to  the  minor  axis  divided  by  sin  e.  Now,  beginning 
with  the  footprint  at  the  point  of  irrpact  for  8 ■ .1,  a set  of  joint 
numbers  bounded  by  the  footprint  ellipse  can  be  determined  for  each  8 
and  listed  on  the  data  code  12  cards.  Footprint  overlap  is  likely  to 
occur,  so  the  same  joint  number  will  then  appear  in  two  or  more  foot- 
print lists.  Nevertheless,  a complete  joint  list  must  be  input  for 
each  footprint  location. 
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The  Loads 
calculates  td 


Generator  uses  data  codes  2,  11  and  12  input  data  and 

from  Equation  (15)  and  F and  F from 

3vg  p 


m v sin  8 
avg  = td 


(19) 


F 


6 


f0  Favg 


(20) 


In  Equation  (19)  the  program  effectively  calculates  sin  fl  from  vector 
operations  involving  b,  n and  d,  so  0 is  not  input.  Having  the  values 
of  Fg  and  the  associated  joint  lists  for  each  increment  the  program 
generates  equilibrium  equations  and  solves  for  a set  of  joint  loads 
which  are  equivalent  to  Fg  for  each  increment.  Incremental  joint 
loads  are  found  from 


V* 


0-1 


(21) 


where  (|)g  are  the  joint  loads,  5(?e  are  the  incremental  loads  and 
3 = 1 ,2,3,...  ,n^s+l.  Finally  the  program  converts  the  incremental 
joint  loads  into  their  components  in  the  joint  X,Y,Z  degrees  of  free- 
dom and  outputs  them  in  the  form  of  a matrix  having  3nj  rows  and  njfi+l 
col  urns.  The  incremental  loads  matrix  is  named  DP0T. 


INITIAL  GENERATOR 

The  input  file  is  Fortran  logical  unit  5 and  is  passed  to  the 
Initial  Generator  from  an  UPDATE  data  file,  as  indicated  in  Figure  10. 
Data  code  decks  required  to  be  present  In  the  input  file  are  shown  in 
Table  4.  A required  data  code  deck  may  consist  of  only  an  end  card 
if  that  type  of  data  is  not  used  In  the  definition  of  a particular 
model , but  the  end  card  must  be  present  as  a minimum.  Optional  data 
codes  present  In  the  Input  will  be  read  over  and  ignored.  The  output 
file  Is  Fortran  logical  unit  1 and  contains  structural  matrices  and 
Information  lists  needed  in  subsequent  steps  of  the  analysis.  Table  6 
shows  the  contents  of  the  output  data  file.  The  output  is  stored  in 
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TABLE  6.  INITIAL  GENERATOR  OUTPUT  FILE 


FORMAT  Tape  Name  Is  MTAPE.l 

NAME 

OF  ARRAY 

DIMENSIONS 

OF  ARRAY 

DESCRIPTION 

UZER0 

nj  x 3 

Matrix  of  joint  coordinates. 

MPT 

n„  x 1 
m 

Table  of  material  properties. 

ECT 

100  x na 
e 

Table  of  element  constants. 

PRUPT 

nu  X nT 

Matrix  of  force  components  in  the  reordered 
unconstrained  global  degrees  of  freedom 
resulting  from  unit  forces  in  the  total 
global  degrees  of  freedom. 

PRUF 

nu  x nPF 

Matrix  of  force  components  In  the  reordered 
unconstrained  global  degrees  of  freedom 
resulting  from  unit  element  forces  in  the 
global  degrees  of  freedom. 

KEL 

nPF  X nPF 

Matrix  of  unassembled  element  stiffnesses 
in  the  element  global  degrees  of  freedom. 

MEL 

nPF  X nPF 

Matrix  of  unassembled  element  masses  in  the 
element  global  degrees  of  freedom. 

KB  EL 

np  x np 

Matrix  of  unassembled  element  stiffnesses 
in  the  element  lumped  degrees  of  freedom. 

CBEL 

Op  x n^ 

Matrix  of  unassembled  element  damping 
sub-matrices  in  the  element  lumped  degrees 
of  freedom. 

FFBAR 

np  x n^ 

Matrix  of  element  forces  In  the  global 
degrees  of  freedom  resulting  from  unit 
element  forces  in  the  element  lumped  degrees 
of  freedom. 

SIGFB 

ns  x n^ 

Matrix  of  average  element  stresses  result- 
ing from  unit  element  forces  in  the  element 
lumped  degrees  of  freedom. 

EPSIG 


n$  x 


Matrix  of  element  strains  resulting  from 
unit  average  element  stresses. 


TABLE  6.  (continued) 


FORMAT  Tape 

s Name  Is  MTAPE.l 

NAME 

OF  ARRAY 

DIMENSIONS 

OF  ARRAY 

DESCRIPTION 

DEBT 

30  x n„ 
e 

Table  of  unassembled  element  deformations 
resulting  from  temperature  gradients. 

EVT 

24  x na 
e 

Table  of  element  variables. 

CONST 

30  x 1 

Table  of  job  constants. 

nj  3 number  of  joints 

n„  3 number  of  words  to  store  material  properties 
m 

ng  3 number  of  elements  excluding  point  mass 

nu  3 number  of  unconstrained  degrees  of  freedom 

3 total  number  of  degrees  of  freedom 
npp  3 number  of  element  forces  including  point  mass 

np  = number  of  element  forces  excluding  point  mass 

np  3 number  of  lumped  element  forces 


n^  3 number  of  element  stresses/strains 


the  form  of  a F0RMAT  matrix  tape  named  MTAPE.l.  It  can  be  passed 
directly  as  an  Input  file  to  the  next  analysis  step  and/or  copied 
to  magnetic  tape  and  saved  for  future  use. 


In  addition  to  generating  all  the  arrays  shown  in  Table  6 and 
creating  the  output  file,  the  Initial  Generator  unconditionally  prints 
out  information  which  may  be  helpful  to  the  user.  Printed  first  is  a 
listing  of  some  of  the  data  code  decks  from  the  input  file..  Included 
are  data  code  2 joint  coordinates,  data  code  4 direction  numbers, 
data  code  5 constraints,  data  code  10  material  properties,  and  data 
codes  20,  30,  40,  50  element  definition  data.  Next  are  two  tables  of 
bookkeeping  Information  and,  finally,  a summary  of  dimensions  and 
input  constants. 

The  first  table  of  bookkeeping  information  defines  the  uncon- 
strained degrees  of  freedom  (UD0F).  There  are  three  translational 
degrees  of  freedom  In  the  global  X,Y,Z  axes  directions  at  each  joint 
in  the  model  and  one  translational  degree  of  freedom  col  inear  with 
each  line  segment  defining  the  edge(s)  of  one  or  more  finite  elements 
in  the  model.  Some  of  the  joint  degrees  of  freedom  may  be  constrained 
as  specified  in  data  code  5 constraints.  Edges  are  not  constrained. 
This  table  lists  all  the  remaining  or  "unconstrained"  joint  ana  edge 
degrees  of  freedom.  Both  joint  UD0F  and  edge  UD0F  are  defined  under 
three  coluim  headings  titled  UD0F,  J1 , and  J2.  The  headings  are 
repeated  five  times  across  the  top  of  a page,  so  that  five  UD0F's 
are  defined  per  line  of  print.  For  joints  the  UD0F  column  contains 
the  degree  of  freedom  Identification  number,  the  J1  column  gives  the 
joint  number,  and  the  J2  column  shows  an  Integer  value  of  T,  2,  or  3 
to  Indicate  the  direction  of  the  degree  of  freedom  is  parallel  to 
the  global  X,  Y,  or  Z axis  of  the  model.  For  example,  the  following 
Information  defines  unconstrained  degree  of  freedom  number  1291  at 
joint  452  in  the  Y-dl recti  on. 


For  edges  the  UD0F  column  contains  the  unconstrained  degree  of  free- 
dom identification  number,  the  J1  and  J2  coluims  give  two  joint  numbers 
defining  a line  segment  (or  edge  of  elements)  which  is  colinear  with 
the  degree  of  freedom  directed  from  J1  toward  J2.  For  example,  the 
following  information  defines  unconstrained  degree  of  freedom  number 
3011  colinear  with  the  line  segment  (edge)  between  joints  164  and  168 
and  directed  from  164  to  168. 

UO0F  J1  J2 

3011  164  -168 

The  negative  sign  attached  to  the  second  joint  number  indicates  that 
the  UD0F  defined  is  an  edge  degree  of  freedom. 

The  second  table  of  bookkeeping  information  gives  a reordered 
numbering  sequence  for  the  unconstrained  degrees  of  freedom.  The 
reordered  sequence  defines  the  row  format  of  matrices  PRUPT  and  PRUF, 
see  Table  6.  In  this  table  the  original  UD0F  identification  numbers 
are  listed  in  columns  titled  "U"  and  the  reordered  numbers  follow 
iirmediately  in  a second  column  titled  "RU".  Ten  pairs  of  columns  are 
printed  across  the  page. 

The  final  summary  print  is  simply  a list  of  the  total  numbers  of 
joints  nj,  constraints  n^,  materials  nM,  bar  elements  nBg,  membrane 
elements  n^,  cell  elements  nCfi,  point  mass  elements  npe>  and  edge 
degrees  of  freedom  n£  together  with  constant  data  including  the  base 
temperature,  mechanical  equivalent  of  heat,  cell  stiffness  suppression 
coefficient  and  bar  compliance  coefficients  (see  data  code  1,  card  2). 
The  user  can  calculate  the  dimensions  defined  at  the  end  of  Table  6 from 
the  summary  data  through  the  use  of  the  following  equations. 


92 


nT  * 3 nj  + nE 

nu  3 nT  - nC 
ne  “ nBe  + nMe  + nCe 
nPF  * 7 nBe  + 16  nMe  + 36  nCe  + 3 nPe 
nF  3 nPF  “ 3 nPe 
"7  3 2 nBe  + 9 nMe  + 30  nCe 
nS  3 nBe  + 3 nMe  + 12  nCe 

MODAL  TRANSFORMATION 

The  extraction  of  a relatively  small  number  of  vibration  modes 
and  using  them  as  a transform  to  reduce  the  problem  size  is  accomplished 
by  means  of  the  FORMAT  Phase  2 program,  see  Reference  1 for  user  docu- 
mentation. The  procedure  Involves  executing  the  three  FORMAT  steps 
shown  in  Figure  10.  The  input  deck  setups  have  been  previously  described, 
see  "FORMAT  Phase  2 Data". 


Y (22) 


A user  who  is  familiar  with  FORMAT  will  recognize  that  he  has 
considerable  freedom  in  setting  up  the  matrix  instruction  sequence 
and  that  he  can  tailor  it  to  suit  his  own  purposes.  For  example,  all 
three  FORMAT  steps  could  be  combined  into  one  in  the  case  of  a small 
finite  element  model  requiring  little  computer  time  and  resources. 
However,  the  three  steps  shown  are  suitable  for  both  large  and  small 
models  and  may  be  executed  as  Individual  computer  runs,  or  optionally 
for  small  models,  linked  by  JCl  and  executed  In  one  computer  run.  If 
the  user  wishes  to  avoid  handling  the  FORMAT  setups  as  card  decks,  he 
may  store  them  In  UPDATE  data  files.  He  may  then  access,  edit  and 
execute  them  as  desired  by  the  use  of  appropriate  JCL. 
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Stiffness  Decomposition 

FORMAT  deck  setup  1 requires  two  Input  files  and  creates  one 
output  file,  reference  Figure  10.  The  first  input  file  is  Fortran 
logical  unit  2 and  is  the  output  file  from  the  Initial  Generator. 

It  is  either  copied  directly  by  F0RMAT  from  the -Initial  Generator  or 
copied  from  magnetic  tape  where  It  was  stored  by  the  Initial  Generator. 

The  second  input  file  is  Fortran  logical  unit  3 and  is  the  output  file 
from  the  Loads  Generator.  It  may  come  directly  from  the  Loads  Generator 
or  it  may  be  copied  from  tape. 

In  the  discussion  of  the  Loads  Generator  output  file  it  was 
pointed  out  that  the  incremental  loads  matrix  DP0T  could  optionally 
be  punched  out  on  cards.  In  that  case  the  user  could  delete  card  4 
(LTAPE.l)  from  the  F0RMAT  deck  setup  1,  and  insert  the  ^MATRIX  card 
followed  by  the  matrix  DP0T  card  data  and  the  "E"  card  all  between  cards 
17  and  18.  Alternately,  instead  of  inputting  the  incremental  loads 
matrix  from  the  Loads  Generator,  the  user  may  elect  to  generate  the 
DP0T  matrix  himself  and  Input  it  on  cards.  Note  that  for  a user 
generated  load  matrix,  the  global  X,  Y,  Z components  of  a load  acting 
on  joint  number  j are  input  in  row  numbers  3j-2,  3j-l  and  3j,  respectively. 
The  row  dimension  of  matrix  DP0T  Is  three  times  the  total  number  of  joints 
and  the  column  dimension  Is  equal  to  the  number  of  incremental  load  changes 
the  user  wishes  to  define.  Finally,  the  user  should  recognize  that  he  can 
use  this  approach  to  combine  bird  impact  incremental  loads  from  the  Loads 
Generator  with  additional  Incremental  loads  which  he  knows  to  be  acting 
on  the  structure  during  impact. 

Note  that  there  Is  a variable  called  3nj  In  F0RMAT  deck  setup  1 
which  must  have  an  Integer  value  assigned  by  the  user.  The  data  card 
Involved  must  be  updated  with  the  correct  value  (3  times  the  total 
number  of  joints  In  the  model)  before  executing  the  job. 

The  output  data  file  Is  Fortran  logical  unit  4 and  contains  the 
following  matrices. 


94 


1)  DP0T  — matrix  of  incremental  loads  representing  the  bird  impact 

2)  MR  — upper  triangular  structural  mass  matrix 

3)  LTl  --  decomposed  stiffness  matrix 

4-)  UR  — static  displacements  in  the  reordered  unconstrained  degrees 
of  freedom  resulting  from  applying  the  first  column  of 
DP0T  as  a set  of  static  loads 

The  output  file  may  be  copied  to  the  appropriate  Fortran  logical 
unit  for  direct  input  to  the  next  FORMAT  step,  and/or  it  may  be  copied 
to  magnetic  tape  and  saved  for  future  use. 

Printout  from  this  step  includes  the  standard  FORMAT  print  and 
matrix  UT,  which  contains  the  static  displacements  of  matrix  UR 
transformed  to  the  more  convenient  total  (T)  degree  of  freedom  format. 
In  the  "T"  format,  global  X,  Y-,  Z displacements  for  joint  number  j 
are  found  in  the  row  numbers  3j-2,  3j-l  and  3j,  respectively.  The 
printout  of  the  UT  matrix  should  be  scanned  to  determine  if  the  model 
is  responding  statically  in  a reasonable  fashion  and  for  any  insights 
it  might  give  as  to  expected  dynamic  response. 


Mode  Extraction  and  Transformation 

F0RMAT  deck  setup  2 requires  two  input  files  and  creates  one  output 
file,  see  Figure  10.  The  first  Input  file  Is  identical  to  that  of 
deck  setup  1.  The  second  Input  file  is  Fortran  logical  unit  3 and  is 
output  from  step  1.  It  Is  either  copied  directly  from  output  unit  4 
in  step  1 to  Input  unit  3 In  this  step,  or  copied  to  unit  3 from 
magnetic  tape  where  It  was  stored  by  step  1. 

The  output  file  is  Fortran  logical  unit  4 and  contains  the  follow- 
ing matrices. 
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1)  TR  — modes  matrix  containing  normalized  displacements  in  the 

reordered  unconstrained  degrees  of  freedom  for  the  first 
n^  natural  vibration  modes 

2)  PBUPJ  — transformed  partition  of  PRUPT  containing  model  displace- 

ments resulting  from  unit  forces  in  the  joint  "T"  degrees  of 
freedom 

3)  PBUF  — transformed  PkUF  matrix  containing  modal  displacements 

resulting  from  unit  element  forces 

4)  DPB0  --  transformed  incremental  loads  matrix 

The  output  file  may  be  copied  to  the  appropriate  Fortran  logical 
unit  for  direct  input  to  the  next  FORMAT  step,  and/or  it  may  be 
copied  to  magnetic  tape  and  saved  for  future  use. 

Note  that  there  is  a variable  called  n^  in  FORMAT  deck  setup  2 
which  must  have  an  integer  value  assigned  by  the  user.  The  data  card 
involved  must  be  updated  with  the  correct  value  before  executing  the 
job.  The  value  n^  is  the  number  of  transformation  modes  desired  by 
the  user.  There  is  no  fixed  rule  for  choosing  the  number  of  modes, 
but  for  practical  purposes  ng  < 50  may  be  considered  as  a guideline. 

In  linear  analysis  of  models  having  5000  degrees  of  freedom,  the  use 
of  30  modes  appears  to  be  adequate.  For  small  models  5 or  10  modes 
may  be  sufficient.  The  core  requirement  can  be  considered  the  same 
as  that  required  for  decomposition,  see  equations  12  and  13. 

In  nonlinear  analysis  limited  experience  Indicates  that  more 
inodes  are  required  and  that  some  of  them  should  be  higher  frequency 
inodes  which  are  not  readily  available.  Suitable  modes  need  to  be 
Identified  and  methods  to  extract  them  need  to  be  developed.  See  "Part  1 
Theory  and  Applications"  (a  separate  volume  of  this  report)  for  further 
discussion. 
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In  addition  to  the  standard  FORMAT  listing  the  printout  from  this 
step  includes  matrices  VAL  and  M0GES.  VAL  is  a column  matrix  contain- 
ing the  reciprocals  of  the  squares  of  the  natural  frequencies  associ- 
ated with  the  vibration  modes  in  matrix  TR.  M0DES  contains  the  normal- 
ized displacement  modes  of  matrix  TR  transformed  to  the  "T”  degree  of 
freedom  format.  To  reduce  the  volume  of  printout  only  values  l .05 
are  printed.  The  frequencies  may  be  of  general  interest  to  the  user. 
The  modes  should  be  plotted  as  they  are  the  deformed  shapes  which  will 
be  combined  to  produce  the  dynamic  response. 


Data  Files  for  Incremental  Solution 

FORMAT  deck  setup  3 requires  two  input  files  and  creates  two  output 
files,  see  Figure  10.  The  first  Input  file  is  Identical  to  that  of 
deck  setups  1 and  2.  The  second  input  file  is  Fortran  logical  unit  3 
and  is  output  from  step  2.  It  is  either  copied  directly  from  output 
wit  4 in  step  2 to  input  unit  3 in  this  step,  or  copied  to  unit  3 
from  magnetic  tape  where  it  was  stored  by  step  2. 

The  output  files  are  Fortran  logical  units  4 and  7.  Unit  4 con- 
tains matrix  PBUF  copied  from  input  file  unit  3 and  matrices  MPT,  UZER0, 
ECT,  MEL,  KBEL , C8EL,  FFBAR,  SIGFB,  EPSIG,  DEBT,  EVT  and  C0NST  copied 
from  input  file  unit  2.  Unit  7 contains  matrices  DPBEf  and  PBUPJ  copied 
from  input  unit  3.  The  output  files  may  be  copied  to  Fortran  logical 
units  20  and  21  for  direct  input  to  the  Incremental  Solution  program, 
and/or  they  may  be  copied  to  magnetic  tape  and  saved  for  future  use. 

The  only  printout  from  this  step  is  the  standard  F0RMAT  print  which 
lists  the  user  input  cards  SF0RMAT  through  ?END  and  prints  the  program 
generated  logic  for  the  solution  of  the  matrix  equations. 
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INCREMENTAL  SOLUTION 

The  Incremental  Solution  program  calculates  the  transient  dynamic 
modal  response  to  the  bird  impact  and  the  resulting  element  forces, 
stresses  and  strains  in  the  finite  element  model.  This  information 
is  calculated  at  each  time  point  and  is  saved  in  the  output  data  file. 

The  program  requires  two  input  data  files  and  creates  one  output  file 
as  shown  in  Figure  10. 

The  input  files  are  Fortran  logical  units  20  and  21.  They  are 
generated  and  output  by  FORMAT  step  3 on  logical  units  4 and  7 and 
are  copied  to  units  20  and  21  for  use  in  the  Incremental  Solution 
either  directly  from  units  4 and  7 or  from  magnetic  tape  where  they 
were  stored  by  FORMAT. 

The  output  files  for  the  Linear  and  Nonlinear  Incremental  Solutions  are 
on  Fortran  logical  units  2 and  23,  respectively.  This  file  contains  the 
information  required  by  the  Postprocessor  program  to  selectively  calculate 
and  print  the  incremental  response  results.  The  output  file  is  generated 
on  magnetic  tape  and  may  be  passed  directly  to  the  Postprocessor  for  initial 
printout  selections.  Also,  it  should  be  saved  for  additional  Postprocessor 
runs  which  the  user  may  want  to  make  after  seeing  the  initial  printout. 

Card  input  data  is  described  under  "Incremental  Response  Data". 

It  consists  of  one  card  of  option  and  dimensional  data  and  one  or 
more  cards  listing  the  time  points  marking  the  end  of  each  increment. 

It  is  important  that  these  time  increments  be  the  same  as  those  used 
in  calculating  the  impact  loads,  see  "Loads  Generator".  The  deck 
setup  for  execution  of  the  Incremental  Solution  code  consists  of  the 
appropriate  JCL  cards  together  with  the  card  input  data  mentioned 
above. 

There  are  two  Incremental  Solution  program  modules  available. 

One  calculates  the  linear  incremental  response;  the  other  calculates 
the  nonlinear  incremental  response.  The  user  may  attach  either  the 
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linear  or  the  nonlinear  program  in  a given  application.  However,  as 
indicated  in  "Mode  Extraction  and  Transformation",  the  use  of  the  non- 
linear program  may  result  In  unrealistic  answers  unless  special  modes 
have  been  used  in  the  modal  transformation.  Further  development 
needs  to  be  done  to  define  and  extract  these  special  modes.  Therefore, 
the  user  should  exercise  caution  in  the  use  of  the  nonlinear  program. 

If  continued  long  enough,  the  dynamic  response  will  eventually  damp 
out  to  some  condition  of  static  equilibrium.  Therefore,  IMPACT  may  be 
used  to  perform  linear  or  nonlinear  static  analyses.  This  use  will  be 
expedited  If  arbitrarily  high  damping  to  stiffness  ratios  are  Introduced 
in  the  material  properties.  Then  a single  relatively  long  time  Increment 
will  suffice  In  the  case  of  a linear  solution.  For  the  nonlinear  solution, 
a number  of  short  time  increments  (as  determined  by  experience)  will  be 
needed  to  reach  static  equilibrium. 

Printout  from  the  Incremental  Solution  program  consists  of  some 
Intermediate  data  generated  during  the  course  of  the  solution  which 
has  been  useful  during  development  of  the  program.  It  Is  In  terms 
of  modal  responses  and  has  limited  physical  significance.  Further 
details  may  be  found  in  Part  3 Programming  Manual. 


POSTPROCESSOR 

The  Postprocessor  program  allows  the  user  to  scan  the  large  amount 
of  output  from  the  Incremental  Solution  and  print  selected  data  In  a 
convenient  format.  It  will  normally  be  run  more  than  once  In  a given 
analysis,  because  after  a first  look  at  a small  amount  of  data  the  user 
will  probably  want  more  detailed  Information  at  critical  locations  In 
the  structure. 

The  Input  file  Is  Fortran  logical  unit  2 and  Is  copied  from  the 
output  tape  generated  by  the  Incremental  Solution.  There  Is  no  out- 
put file.  Printout  Is  controlled  by  card  Input  data  as  described 
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under  "Postprocessor  Data".  The  first  card  allows  the  user  to  specify 
whether  or  not  he  wants  to  exercise  selectivity  with  regard  to  time 
Increments,  joints,  bars,  membranes,  and/or  cells.  It  also  provides 
for  the  input  of  a displacement  scale  factor.  Cards  2 through  6 are 
used  to  indicate  the  type  of  data  to  be  printed,  such  as  displacements, 
forces,  stresses,  etc.  Cards  7 through  10  provide  title  information 
which  will  be  printed  at  the  top  of  each  page  of  printout.  Finally, 
cards  20,  30,  40,  50,  60  contain  user  selected  time  increment  numbers, 
joint  numbers,  bar  numbers,  membrane  numbers,  and  cell  numbers, 
respectively,  consistent  with  the  selectivity  options  Indicated  In 
card  1. 

The  deck  setup  for  execution  of  the  Postprocessor  program  consists 
of  the  appropriate  JCL  cards  together  with  the  card  input  data  discussed 
above.  The  following  pages  show  samples  of  the  printout  formats  from 
the  Postprocessor.  Note  that  the  results  printed  are  at  the  end  of  a 
long  (10  sec)  time  Increment,  and  that  the  velocities  and  accelerations 
have  damped  out  so  that  zero  values  are  printed  out.  The  results  are 
applicable  to  the  finite  element  model  shown  in  Figure  14. 
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Sample  printout  of  joint  displacements  and  coordinates 
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Sample  printout  of  bar  element  forces,  stresses  and  strains 
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Sample  printout  of  membrane  element  forces,  stresses  and  strains 
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Sample  printout  of  cell  element  stresses  and  strains 


SECTION  IV 

PROGRAM  DIAGNOSTICS  AND  ERROR  MESSAGES 

This  section  contains  lists  of  program  diagnostics  and  error  messages 
for  the  IMPACT  Programs.  In  general,  the  messages  are  designed  to  point 
to  discrepancies  in  the  input  data  and  are  self-explanatory.  Most  repre- 
sent fatal  errors,  although  the  program  will  continue  to  check  as  much 
of  the  Input  data  as  possible  before  the  termination  of  the  run.  Some 
of  the  messages  are  informational  or  call  the  user's  attention  to  a 
possible  discrepancy,  but  they  do  not  terminate  execution.  In  some 
cases  a message  may  indicate  Insufficient  space  allocation.  Then  the 
user  may  have  to  refer  to  Part  3 Programming  Manual  and/or  seek  program 
maintenance  help. 

LAMINATE  GENERATOR 

The  following  is  a list  of  diagnostic  messages  printed  by  the  Laminate 
Generator  program.  They  are  fatal  errors  unless  noted  otherwise. 

1.  NO  LAYERS  (DATA  COOE  7)  ENCOUNTERED  IN  INPUT.  JOB  TERMINATED. 

2.  NO  CELLS  (DATA  COOE  40)  ENCOUNTERED  IN  INPUT.  JOB  TERMINATED. 

3.  NO  JOINTS  (DATA  COOE  2)  ENCOUNTERED  IN  INPUT.  JOB  TERMINATED. 

4.  JOINT  NO.  XXXX  LISTED  IN  DATA  CODE  8 IS  ABOVE  A USER  DEFINED  JOINT  - 
USER  DEFINED  JOINT  WILL  REMAIN  UNCHANGED,  (non-fatal) 

5.  CARO  ENCOUNTERED  IN  DATA  CODE  8 DUPLICATES  EXISTING  LAMINATE  OATA. 
(non-fatal) 

6.  JOINT  NO.  XXXX  LISTED  IN  DATA  CODE  8 IS  NOT  FOUND  IN  LIST  OF  QUAD 
CORNER  JOINTS  - BAD  JOINT  OATA  MAY  RESULT,  (non-fatal) 

7.  NO  REVISED  THICKNESSES  (DATA  CODE  8)  ENCOUNTERED  IN  INPUT,  (non-fatal) 
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3. 


GAP  IN  SEQUENCE  (non-fatal) 


9.  ATTEMPT  TO  UNITIZE  ZERO  VECTOR,  (non-fatal) 

LOADS  GENERATOR 

The  following  is  a list  of  diagnostic  messages  printed  by  the  Loads 
Generator  Program.  They  are  fatal  unless  noted  otherwise. 

1.  A 29999  CARD  WAS  NOT  FOUND. 

2.  NO  JOINT  COORDINATE  DATA  (DATA  CODE  2)  ENCOUNTERED  IN  INPUT. 

3.  NO  CONSTANTS  (DATA  CODE  11)  ENCOUNTERED  IN  INPUT. 

4.  INCREMENT  (BETA)  NUMBER  XXXX  NOT  FOUND  IN  DATA  CODE  11  DATA. 

5.  NO  LOADED  JOINT  NUMBERS  (DATA  CODE  12)  ENCOUNTERED  IN  INPUT. 

6.  JOINTS  FOR  INCREMENT  (BETA)  NUMBER  XXXX  NOT  FOUND  IN  DATA  CODE  12  DATA. 

7.  ATTEMPT  TO  UNITIZE  ZERO  VECTOR  (non-fatal) 

INITIAL  GENERATOR 

The  following  is  a list  of  diagnostic  messages  printed  by  the  Initial 
Generator  Program.  They  are  grouped  under  subroutine  names  which  are 
not  described  here.  Subroutine  descriptions  may  be  found  in  Part  3, 
Programming  Manual. 

BIR0G1 

1.  DATA  CODE  1 INPUT  IN  ERROR. 

2.  (BIRDG1 ) TBASE  » 0 DATA  INPUT  ERROR.  JOB  ABORTED.  DATA  SCAN  CONTINUES. 
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ACCPRP 


1.  (ACCPRP)  MATERIAL  REFERENCE  NUMBER  OUT  OF  RANGE. 

REQUESTED  MATERIAL  NUMBER  NO.  NNNNN 
MAXIMUM  MATERIALS  NNNNN 

2.  FATAL  ERROR  CONDITION  EXISTS  OUTSIDE  OF  THE  MATERIAL  PROPERTIES  ACCESS 
MODULE , NO  ACCESS  ATTEMPTED. 

CNSTRN 

1.  ALL  DIRECTION  COMPONENTS  ARE  ZERO. 

DATA  CODE  - NN 

CONSTRAINT  NO. * NNNNN 
JOINT  M * NNNNN 

JOINT  N = NNNNN 

DIRECTION  NO.  ■ NNNNN 
DIRECTIONS 
NNNNNNN ,NNN 

2.  INVALID  DATA  IN  THE  DIRECTION  TABLE  ASSOCIATED  WITH  DIREC  NO.  N 

3.  INVALID  DIRECTION  NUMBER.  MAX  DIREC.  NO.»NNNNN 

4.  INVALID  JOINT  NUMBER.. 

5.  CONSTANT  TABLE  DATA  OUT  OF  SEQUENCE,  CONSTRAINT  NUMBER  EXPECTED  NNNNN 

6.  INVALID  DATA  CODE. 

7.  THE  NUMBER  OF  CONSTRAINTS  EXCEEDS  THE  PROGRAM  MAXIMUM  OF  NNNNN 

8.  CONSTRAINT  DATA  INPUT  MISSING. 

9.  THE  NUMBER  OF  OBLIQUE  CONSTRAINTS  EXCEEDS  THE  MAXIMUM  OF  NNNNN 


COCALC 


1 .  ERROR  DETECTED  IN  SUBROUTINE  COCALC 

MATRIX  ASC  IS  EITHER  SINGULAR  OR  NEAR  SO. 

DRCNUM 

1.  ALL  DIRECTION  NUMBERS  ARE  ZERO.  A UNIT  VECTOR  PARALLEL  TO  THE  X AXIS 
IS  ENTERED. 

2.  JOINT  NUMBER  OUT  OF  RANGE. 

MAXIMUM  JOINT  * NNNNN 

JOINT  M - NNNNN 

JOINT  N * NNNNN 

DIRECTION  NUMBER  NNNNN 

FIXUP  TAKEN,  BAD  DATA  FLAGGED  IN  DIREC  TBL. 

3.  DIRECTION  COSINE  RECORDS  NOT  IN  SEQUENCE. 

NNNNN  RECORDS  READ 

ERRONEOUS  RECORD 

ID  NDIREC  JNTM  JNTN  IBLNK  DIREC 

4.  INVALID  DATA  CODE  IN  DIRECTION  NUMBER  INPUT 
ERRONEOUS  RECORD 

ID  NDIREC  JNTM  JNTN  IBLNK  DIREC 

5.  THE  NUMBER  OF  DIRECTION  NUMBERS  EXCEEDS  THE  PROGRAM  MAXIMUM  OF  50. 
EDGDOF 

1.  ***  SUBROUTINE  EDGDOF,  FATAL  ERROR  PRIOR  TO  ENTRY 

EDGDOF  WILL  NOT  BE  EXECUTED. 

2.  ***  ERROR  DETECTED  IN  SUBROUTINE  EDGDOF 

NNNNN  JOINT  PAIRS  WRITTEN  TO  TAPE  IN  SU8R  EDGES 
NNNNN  JOINT  PAIRS  READ  FROM  TAPE  IN  EDGDOF 
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EDGES 

1.  ***  ERROR  IN  SUBROUTINE  EDGES 
INVALID  DATA  CODE  NNNNNNNN 


ELEMNT 


1.  (ELEftrr)  ELEMENT  NO.  OUT  OF  SEQUENCE. 

PREVIOUS  SEQUENCE  NO.  » NNN 
CURRENT  CARD  SEQUENCE  NO.  - NNN 


2. 


(ELEMNT) 


ELEMENT  TYPE  CODE  OUT 
CARD  MISSING, 

PREVIOUS  TYPE  CODE  • 
CURRENT  CARD  TYPE  - 
ELEMENT  NUMBER 


OF  SEQUENCE  OR  9999  SEPARATOR 

NNN 

NNN 

NNNNN 


3.  (ELEWT)  A ■ 0 FOR  BAR  ELEMENT  NNNNN 


4.  (ELEMNT)  T - 0 FOR  MEMBRANE  ELEMENT  NNNNN 


5.  (ELEMNT)  MR  » 0 FOR  ELEMENT  NNNNN  CODE  • NNN 


6.  (ELEMNT)  JOINT  NUMBER(S)  FOR  ELEMENT  NNNNN  IS  (ARE)  GREATER  THAN  MAX 

JOINT  INPUT  ■ NNNNN 

JOINT  NUMBERS  INPUT  ARE  PO,  QO,  RO,  SO,  PI,  Q1 , R1 , SI 

7.  (ELEMNT)  JOINT  NUMBER(S)  FOR  ELEMENT  NNNNN  IS  (ARE)  EQUAL  ZERO, 

JOINT  NUMBERS  INPUT  ARE  PO,  QO,  RO,  SO,  PI,  Q1 , R1 , SI 


JTEMP 

1.  (JTEMP)  DATA  CODE  ERROR  READING  TEMPERATURE  DATA. 

ICODE  ■ NN  J1.J2  ■ NN  NN  TEMP  ■ .NNNNNN  ENN 
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2.  ( JTEMP ) JOINT  NO.  NOT  INPUT  FOR  TEMPERATURE  DATA. 

I CODE  * NN  J1,J2  * NN  NN  TEMP  • . NNNNNN  ENN 

JZERO 

1.  (JZERO)  DATA  CODE  ERROR  READING  COORDINATE  DATA. 

I CODE  =>  NN  J * NN  X,Y,Z  » .NNNNN  ENN  .NNNNN  ENN  . NNNNN  ENN 

2.  (JZERO)'  JOINT  NO.  NOT  INPUT  FOR  COORDINATE  DATA. 

I CODE  = NN  J * NN  X,Y,Z  * .NNNNN  ENN  .NNNNN  ENN  .NNNNN  ENN 

3.  (JZERO)  JOINT  NO.  OUT  OF  SEQUENCE. 

I CODE  =»  NN  J * NN  X,Y,Z  * .NNNNN  ENN  .NNNNN  ENN  .NNNNN  ENN 

LPCKC 

1.  LPCKC-SID  ***  KBAR  = D INVERSE  ERROR. 

CELL  DEFINED  BY  JOINTS  NNNN  NNNN  NNNN  NNNN  NNNN  NNMN  NNNN 
NNNN  HAS  GENERATED  A SINGULAR  MATRIX  AND  CANNOT  BE  INVERTED. 

2.  LPCKC-SID  ***  KBAR  = D MATRIX  HAS  ONLY  ABOUT  NNNN . N SIGNIFICANT 
DIGITS. 

CELL  DEFINED  BY  JOINTS  NNNN  NNNN  NNNN  NNNN  NNNN  NNNN  NNNN 
NNNN 

LPMKC 

1.  LPMKC-SID  ***  KBAR  - D INVERSE  ERROR. 

MEMBRANE  DEFINED  BY  JOINTS  NNNN  NNNN  NNNN  NNNN 
HAS  GENERATED  A SINGULAR  MATRIX  AND  CANNOT  BE  INVERTED. 

2.  LPMKC-SID  ***  KBAR  - D MATRIX  HAS  ONLY  ABOUT  NNNN . N SIGNIFICANT 
DIGITS. 

MEMBRANE  DEFINED  BY  JOINTS  NNNN  NNNN  NNNN  NNNN 
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MATDES 


1 . MATERIAL  REFERENCE  NUMBER  NN  BEING  PROCESSED 
AND  MATERIAL  REFERENCE  NO.  NN  WAS  ENCOUNTERED. 

THE  ERRONEOUS  RECORD 

IDCODE  NTRY  NPC  MAP  REC  COHT 

***  ERROR  DETECTED  IN  SUBROUTINE  MATDES 

2.  INVALID  IDCODE  NN  ENCOUNTERED  WHILE  PROCESSING  MATERIAL 
PROPERTIES  DATA 

THE  ERRONEOUS  RECORD  . 

IDCODE  NTRY  NPC  MAP  REC  COHT 

***  ERROR  DETECTED  IN  SUBROUTINE  MATDES 

3.  INSUFFICIENT  SPACE  ALLOCATED  FOR  MATERIAL  DESCRIPTION 
STATEMENTS,  MODIFY  VARIABLE  MXDES  IN  SUBR.  MTU  10 D. 

THE  ERRONEOUS  RECORD  . 

IDCODE  NTRY  NPC  MAP  REC  CONT 
***  ERROR  DETECTED  IN  SUBROUTINE  MATDES 

MTLMOD 

1 . ***  ERROR  DETECTED  IN  SUBROUTINE  MTLMOD 

NO  MATERIAL  DESCRIPTION  INPUT  FOR  REFERENCE  NUMBER  NWNNN 

2.  ***  ERROR  DETECTED  IN  SUBROUTINE  MTLMOO 
MATERIAL  REFERENCE  NO.  OUT  OF  RANGE 

SPACE  ALLOCATED  FOR  NNNNN  ERRONEOUS  MTRL . REF.  NO.  NNNNN 
IF  MODIFICATION  OF  THE  PROGRAM  IS  NECESSARY  THEN 
MODIFY  VARIABLE  MXMTRL  IN  SUBR.  MTLMOO, 

MODIFY  SUBSCRIPTS  IN  VARIABLES  LNGTHS , 

AND  MPTR  IN  LABELED  COMMON  MTRL. 

3.  ***  ERROR  OETECTED  IN  SUBROUTINE  MTLMOO 
INVALID  DATA  CODE  NN  ENCOUNTERED 
IDCODE  NTRY  NPC  MAP  REC  CONT 
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4.  ***  ERROR  DETECTED  IN  SUBROUTINE  HTLMOD 
MATERIAL  REFERENCE  NO.  OUT  OF  SEQUENCE  OR  MISSING 
MATERIALF  REFERENCE  NO.  LAST  ENCOUNTERED  NNNNN 
MATERIAL  REFERENCE  NO.  JUST  ENCOUNTERED  NNNNN 

5.  ***  ERROR  OETECTED  IN  SUBROUTINE  MTLMOD 
0VER3 

1.  (0VER3)  ERROR  IN  CONSTRAINT  TABLE  GENERATION 

MORE  THAN  3 ENTRIES  FOR  JOINT  NNNNN 
CONSTRAINT  NO.  JOINT  NO.  DIRECTION  VALUE 
NNNNN  NNNNN'  NNNNN 

PASSM2 

1.  ERROR  (PASSM2)  DUPLICATE  CONSTRAINTS  AT  JOINT  NNNNN 
PROP 

1.  PROCESSING  MATERIAL  REF.  NO.  NNNNN  PROPERTY  NO.  NNNNN 
SPACE  ALLOCATED  FOR  NNNNN  PROPERTIES  MAXIMUM. 

MODIFY  VARIABLE  MXPROP  IN  SUBR.  MTLMOD, 

MODIFY  SUBSCRIPTS  IN  VARIABLES  LNGTHS , MPTR,  AND 
PRPVAL  IN  LABELED  COMION  MTRL. 

MOD.  SUBSCRIP  IN  VAR.  COEF  AND  VAR  MXPROP  IN  SUBROUTINE  ACCPRP 
***  ERROR  OETECTED  IN  SUBROUTINE  PROP 

2.  PROCESSING  MATERIAL  REF.  NO.  NNNNN  PROPERTY  NO.  NNNNN 
THE  NUMBER  OF  COEF.S  READ  NNNNN  DOES  NOT  EQUAL  THE  NUMBER 
SPECIFIED  TO  BE  INPUT  NNNNN 

3.  PROCESSING  MATERIAL  REF.  NO.  NNNNN  PROPERTY  NO.  NNNNN 
THE  NUMBER  OF  COEF.S  SPECIFIED  TO  BE  CALCULATED  NNNNN  IS 
GREATER  THAN  THE  NUMBER  OF  VALUES  INPUT  NNNNN 
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4.  PROCESSING  MATERIAL  REF.  NO.  NNNNN  PROPERTY  NO.  NNNNN 
THE  NUMBER  OF  COEF.S  NNNNN  WILL  OVERFLOW  THE  SPACE  NNNNN 
ALLOCATED  FOR  ALL  COEFFICIENTS. 

MODIFIED  VARIABLE  MXCOFS  IN  SUB.  MTLMOD , 

MODIFY  SUBSCRIPTS  IN  VARIA3LE  COEFS  IN  LABELED  COMMON  MTRL. 

5.  PROCESSING  MATERIAL  REF.  NO.  NNNNN  PROPERTY  NO.  NNNNN 

MORE  PROPERTY  VALUES  TO  BE  INPUT  THAN  SPACE  AVAILABLE  NNNNN 
MODIFY  VARIABLE  MXPVL  IN  SUB.  MTLMOD, 

MODIFY  SUBSCRIPTS  IN  ARRAYS  TVPA,  TVPAT,  AND  VPT. 

6.  PROCESSING  MATERIAL  REF. NO.  NNNNN  PROPERTY  NO.  NNNNN 

MORE  COEFFICIENTS  TO  BE  INPUT  NNNNN  THAN  SPACE  AVAILABLE  NNNNN 
MODIFY  VARIABLE  MXPCO  IN  SUB.  MTLMOD, 

MODIFY  SUBSCRIPTS  IN  ARRAYS  TVPA,  TVPAT,  AP,  RTMP, 

ITMP,  ASC,  AND  BSC  IN  LABELED  COMMON  COFCLC. 

7.  PROCESSING  MATERIAL  REF. NO.  NNNNN  PROPERTY  NO.  NNNNN 
MATERIAL  REFERENCE  NUMBER  WAS  EXPECTED. 

8.  DATA  CODE  NNNNN  WAS  ENCOUNTERED  INSTEAD  OF  10. 

ERRONEOUS  RECORD 

IDCODE  NTRY  NPC  MAP  REC1  REC2  REC3  REC4  REC5  REC6 
RAMOSG 

1.  (RAMOSG)  R * 1.0 

MTLERR  EA/E  SR  R M K «1 

.NNNN  ENN  .NNNN  ENN  .NNNN  ENN  . NNNN  ENN  .NMNN  ENN 
MR  E EA  SIGA  SIGR 

NN  .NNNN  ENN  .NNNN  ENN  .NNNN  ENN  .NNNN  ENN 

2.  (RAMOSG)  ZERO  ARGUMENT  CALCULATED  FOR  NATURAL  LOG  FUNCTION. 

MTLERR  » HNN 

MR  E EA  SIGA  SIGR 

NN  .NNNN  ENN  .NNNN  ENN  .NNNN  ENN  .NNNN  ENN 

3.  (RAMOSG)  RAMBERT-OSGOOO  VALUES  DID  NOT  CONVERGE  WITHIN  50 

1 ITERATIONS,  MTLERR  - NN. 

i MR  E EA  SIGA  SIGR  J 

NN  .NNNN  ENN  .NNNN  ENN  .NNNN  ENN  .NNNN  ENN  1 


J 
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SKDATA 

1.  (SKDATA)  INPUT  DATA  OUT  OF  SEQUENCE,  ATTEMPT  TO  SKIP  CODE  NN 
FAILED  WHEN  THE  FOLLOWING  CARD  WAS  ENCOUNTERED. 

ID9999A  A A A A A 

N N Y Y Y AAAAAAAAAA A AAAAA AAA AAA AAAAAAAAAA AA A A AA AA AA A AAA A A 

FORMAT  PHASE  2 

Diagnostic  messages  printed  by  the  FORMAT  program  are  not  listed  here 
but  are  available  in  Reference  1. 


INCREMENTAL  SOLUTION 

The  following  is  a list  of  diagnostic  messages  printed  by  the  Incremental 
Solution  program: 

1.  TIME  NOT  SPECIFIED  IN  ASCENDING  ORDER.  STOP. 

2.  INSUFFICIENT  WORK  SPACE  3Y  XXXXX  WORDS.  STOP. 

3.  UNCORRECTABLE  SINGULAR  EIGENVECTORS.  HALT. 

POSTPROCESSOR 

When  an  error  is  detected  by  the  Postprocessor  program,  the  following 
message  is  printed  at  the  top  of  the  last  page  of  output. 

*****  ERROR  nnnn  DETECTED  ***** 

The  error  number  nnnn  is  an  integer  having  the  values  and  meanings  described 


nnnn 


Description 


* 


1-6  Input  card  nnnn  is  out  of  sequence 


7 

21 

22 

30 

31 

32 


40 


41 


42 

50,51 ,52 
60,61 ,62 
100 


One  of  the  input  title  cards  (7  through  10)  is  out  of 
sequence 

Increment  selectivity  option  was  specified  but  selection 
table  (card  20)  is  not  present 

Null  Increment  selection  table 

Inconsistent  options  for  joint  data  between  selectivity 
flag  and  joint  subset  selection  flags 

Joint  selectivity  option  was  specified  but  selection  table 
(card  30)  is  not  present 

Null  joint  selection  table 

Inconsistent  options  for  bar  data  between  bar  selectivity 
flag,  bar  subset  selection  flags  and/or  presence  of  bar 
elements  in  model 

Bar  selectivity  option  was  specified  but  selection  table 
(card  40)  is  not  present 

Null  bar  selection  table 

Same  as  40,  41,  42  except  substitute  "membrane"  for  "bar" 
Same  as  40,  41,  42  except  substitute  "cell"  for  "bar" 
Insufficient  core  available 


When  the  error  number  nnnn  is  equal  to  100,  the  following  additional 
message  is  printed  out. 

INSUFFICIENT  CORE  XXXX  VS  YYYY 

The  first  Integer  value  XXXX  Is  equal  to  the  core  available  and  the  second 
value  YYY  Is  the  core  required. 


I 
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SECTION  V 


ILLUSTRATIVE  EXAMPLES 


Figure  14  shows  a simple  finite  element  model  composed  of  two  bar  elements, 
two  membrane  elements  and  two  cell  elements.  Its  orientation  with  respect  to 
the  global  coordinate  axes  is  shown,  joint  and  element  numbering  sequences 
are  given,  and  a bird  Impact  loading  condition  is  defined.  Table  7 illustrates 
the  use  of  the  key-punch  data  forms  to  list  the  model  definition  data  in  the 
proper  card  formats.  Note  that  data  codes  2 and  40  are  not  complete  and  that 
data  codes  7 and  3 have  been  provided  for  a laminate  generator  run  which  will 
calculate  the  missing  data. 


The  data  code  3 card  specifies  a constant  temperature  at  all  joints,  which 
is  also  equal  to  the  base  temperature.  Therefore,  there  will  be  no  thermal 
deformations  in  the  model.  The  data  code  4 cards  define  three  direction 
numbers  parallel  to  the  X,  Y,  Z coordinate  axes.  Data  code  5 cards  impose 
the  constraints  indicated  in  Figure  14.  Notice  that  the  three  available  options 
are  all  exercised  in  defining  the  constraint  directions.  In  this  case  the 
most  convenient  option  is  probably  the  reference  to  data  code  4 direction 
numbers . 


The  model  in  Figure  14  has  only  one  layer  of  cell  elements,  but  the 
Laminate  Generator  can  still  be  used  to  calculate  lower  surface  joint  coor- 
dinates and  to  complete  the  single  layer  cell  definitions.  In  the  data  code 
7 cards  a laminate  number  of  100  Is  specified,  and  the  surface  offset  field 
Is  blank  (tQ  ■ 0).  Therefore,  the  program  will  use  the  surface  joint  coor- 

inates  as  given  without  offset  normal  to  the  surface.  Only  one  layer 
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PRECEDING  PAGE  BUMC-NOT  fILMKQ 


' 3 Constraint 
B1  ,B2  * Bars 
Ml , M2  * Membranes 
Cl,  C2  * Cells 
sin  9 * .6 
Bird  mass  = .004,  Velocity  = 5000,  L = 15 
Length  of  time  Increments  3 .001  sec 
A-j , Ag,  Aj  ■ Travelling  footprint  c.g.'s 


Figure  14.  A Slmole  Finite  Element  Model  for  Purposes  of  Illustration. 


definition  card  (L  + 1 • 101)  is  needed  in  this  case,  but  if  the  laminate 
had  more  than  one  layer  there  would  be  additional  layer  definition  cards  102, 

103,  etc.  The  data  code  8 cards  are  not  really  necessary  since  no  change  in 
layer  thicknesses  is  being  made,  but  they  illustrate  the  procedure.  If  the 
laminate  had  more  layers,  fields  t£,  t-j,  etc.  would  also  contain  overriding 
thickness  values. 

The  data  code  10  cards  illustrate  the  input  of  room  temperature  properties 
for  aluminun  based  on  equations  15  through  18)  and  for  polycarbonate  (based 
on  preliminary  data  for  Reference  3).  Temperature  dependence  for  any  property 
would  be  input  by  specifying  the  desired  number  of  polynomial  coefficients  in 
field  NC  and  giving  temperature,  property  value  pairs  in  fields  T^ , V^. 

Data  codes  11  and  12  are  provided  so  that  a Loads  Generator  run  can  be 

executed  to  calculate  the  bird  impact  loads.  In  this  case  a travelling  foot- 

print loading  has  been  specified  with  a rectangular  force-time  distribution, 
see  data  code  11,  cards  6.  The  user  should  recognize  that  he  may  also  specify 
a stationary  footprint  loading  by  setting  D * 0 for  all  P.  Likewise,  the 
force-time  distribution  may  be  controlled  by  means  of  the  load  distribution 
factors  f. 

Element  definitions  on  data  code  20,  30,  40  and  50  cards  show  complete  bar 
and  membrane  definitions  but  only  partial  cell  data.  The  Laminate  Generator 
will  produce  the  data  for  fields  p-j , q-j , r-j,  s.]  and  M.  If  the  stress  orientation 
angle  £pq  Is  non-zero  for  any  surface  cell,  It  should  be  entered  by  the  user, 

since  the  Laminate  Generator  does  not  calculate  this  angle.  However,  it  does 

duplicate  the  surface  cell  angle  input  by  the  user  for  all  other  cells  in  a 
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stack  in  the  case  of  a multi-layered  laminate. 


Table  8 illustrates  FORMAT  deck  setup  1 card  data  entered  on  a key-punch 
form.  Deck  setups  2 and  3 are  generated  in  a similar  manner. 

Table  9 shows  the  data  cards  for  the  Incremental  Solution  and  Postprocessor 
programs.  The  impact  loading  will  be  acting  on  and  travelling  along  the  sur- 
face of  the  structure  for  the  first  three  increments  only,  but  inertia  forces 
could  cause  more  critical  response  after  impact  than  during  impact.  Therefore, 
ten  time  increments  are  specified  for  the  incremental  solution.  The  Post- 
processor cards  specify  that  results  will  be  printed  for  time  increments  1 
tnrough  6.  Tne  results  printed  are  to  include  the  following: 

1.  Joint  coordinates,  displacements,  velocities  and  accelerations 
for  joints  1 through  12. 

2.  Modal  displacements,  velocities  and  accelerations. 

3.  Lumped  element  forces,  average  stresses,  average  strains  and 
equivalent  stresses  for  all  bar,  membrane  and  cell  elements. 

For  larger  models  than  that  shown  in  Figure  14,  the  user  will  normally 
make  a first  Postprocessor  run  to  print  a more  limited  selection  of  results. 
Then,  after  Inspection  of  this  select  Information,  he  will  know  better  what 
Information  to  print  out  in  subsequent  Postprocessor  runs  for  detailed 
evaluation  of  results. 


i 
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TABLE  7.  MODEL  DEFINITION  DATA  LISTED  ON  KEY-PUNCH  FORMS 


TABLE  7.  (Continued) 
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TABLE  7.  (Continued) 
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TABLE  7.  (Continued) 
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TABLE  7.  (Continued) 
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TABLE  7.  (Continued) 
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TABLE  9.  INCREMENTAL  SOLUTION  AND  POSTPROCESSOR  DATA  CARDS 


APPENDIX  A 

SECTIONS  II  AND  III  OF 
AFFDL-TR-77-99  PART  3 


This  appendix  contains  direct  duplications  Sections  II  and  III  of 
the  Part  3 Programming  Manual.  These  sections  are  presented  here  for  the 

I purpose  of  providing  the  user  with  example  deck  setuDS,  including  JCL,  for 

standard  operations  at  the  Wright-Patterson  Air  Force  Base  CDC  Computing 
Facility.  At  other  installations  this  information  may  not  be  completely 
applicable,  but  additional  information  concerning  imolementation  and 
execution  of  the  IMPACT  Programs  is  available  in  Part  3. 
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SECTION  II 
SYSTEM  OVERVIEW 

The  user  Initiates  an  analysis  by  preparing  tabular  data  which  des- 
cribes the  Idealized  model  of  the  windshield  system  under  Investigation. 
The  model  may  Include  support  structure  for  the  transparency  which  may 
Itself  consist  of  a number  of  laminates.  The  tabular  data  prepared  by 
the  user,  consisting  of  geometry,  constraint,  material  property,  and 
element  definition  data.  Is  described  In  detail  In  Part  2. 

PROGRAM  FUNCTIONS 

The  Laminate  Generator  provides  a preprocessing  function  with  respect 
to  multi-layered  transparencies*  Optionally,  the  user  may  prepare  the 
model  definition  data  consistent  with  the  requirements  of  the  Laminate 
Generator.  The  program  will  then  augment  the  Input  model  definition  data 
with  the  joint  coordinates  and  element  definitions  of  Interior  joints  and 
elements  within  a laminate.  In  addition,  surface  normal  vectors  repres- 
enting a unit  pressure  distribution  over  the  exterior  surface  of  the 
transparency  are  generated.  When  applicable,  this  program  reduces  the 
amount  of  Input  data  prepared  by  the  user  and  provides  vector  data  for 
use  In  assembling  pressure  loading  conditions  In  subsequent  analysis 
steps . 

The  Initial  Generator  serves  the  prime  function  of  transforming  the 
complete  model  definition  data  In  tabular  form  Into  matrix  data  consist- 
ent with  the  requirements  of  the  subsequent  solution  process.  Using  the 
model  definition  data,  the  unassembled  element  stiffness,  damping,  and 
mass  matrices  are  generated  as  are  other  matrices  of  geometric  trans- 
formations, thermal  effects,  etc.  The  data  generated  here  Is  Independ- 
ent of  mechanical  loading  on  the  structure  and  satisfies  the  require- 
ments of  either  a linear  and  non-linear  analysis  whichever  might  follow. 

The  Loads  Generator  provides  the  means  of  computing  discrete  joint 
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loads  on  the  Idealized  model  as  a function  of  time  and  bird  length,  mass, 
velocity,  and  Impact  direction.  Optionally,  the  user  may  prepare  the 
model  definition  data  to  Include  the  description  of  the  Impact  loading 
consistent  with  the  requirements  of  the  Loads  Generator.  The  program 
will  then  generate  a matrix  of  Incremental  joint  loads  at  the  times 
specified  in  either  card  image  or  binary  form. 

The  FORMAT  Phase  II  program  is  used  to  assemble  the  structural  mass 
and  stiffness  matrices,  decompose  the  structural  stiffness  matrix,  ex- 
tract vibration  modes,  perform  the  modal  transformation,  and  produce 
properly  ordered  binary  output  data  files  for  Input  to  the  Incremental 
solution  step.  Input  to  this  processing  step  are  the  binary  file  of 
matrices  from  the  Initial  Generator  and,  optionally,  the  loads  matrix 
from  the  Loads  Generator  in  either  card  Image  form  or  as  a binary  matrix 
In  a separate  file.  The  required  tasks  in  this  phase  of  the  solution 
process  may  be  accomplished  in  one  or  more  executions  of  the  FORMAT 
program. 

The  Linear  and  Non-Linear  Incremental  Solution  programs  solve  the 
equations  of  motion  Incrementally  for  the  modal  response,  structural 
displacements,  and  element  forces,  stresses,  and  strains.  The  results 
from  each  Increment  are  output  to  a binary  file  for  subsequent  process- 
ing. The  format  and  content  of  the  output  file  is  identical  for  either 
a linear  or  non-linear  solution. 

The  Postprocessor  accepts  the  output  file  from  either  a linear  or 
non-linear  solution  and  selectively  prints  the  results.  Selection  may 
be  specified  for  time  Increments,  joints,  modes,  and  elements  In  any 
logical  combination. 

DATA  FILES 

Figure  A.l  shows  the  flow  of  data  through  a complete  solution. 

Each  data  file  Is  identified  for  subsequent  reference  in'  this  discussion. 
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Figure  A.l.  Data  Flow  Through  Complete  Solution 
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e.g..  Cl,  M3,  etc.  The  FORTRAN  logical  unit  numbers  of  the  Input  and 
output  files  for  each  program  are  shown  In  parenthesis.  Not  shown  In  the 
figure  Is  system  print  file  6 which  Is  used  for  printed  output  by  each 
program  nor  Is  the  use  of  the  CDC  utility  program  UPDATE  which  Is  dis- 
cussed later  in  this  section. 

File  Formats 

There  are  only  two  file  structures  used  for  all  data  files  shown  in 
Figure  A.l.  Those  files  Identified  with  the  letter  "C”  are  sequential 
files  of  card  Images  (formatted).  Each  logical  record  in  these  files 
corresponds  to  a single  card  and  contains  a string  of  80  characters. 

The  formats  associated  with  each  card  are  given  in  Part  2. 

Those  files  identified  with  the  letter  "M"  are  sequential  files  of 
matrix  data  (binary,  unformatted).  The  structure  of  these  files  corres- 
ponds to  that  of  FORMAT  matrix  data  (Reference  1).  For  purposes  of  con- 
tinuity, a description  of  this  file  structure  is  repeated  here. 

Each  logical  record  in  these  files  is  of  the  form 
IC01,  K0DE,  NUM,  (A(I ) .1-1 ,NUM) 

where  IC0L  is  an  integer  identifier  for  the  record,  K0DE  is  an  integer 
flag  indicating  compressed  or  expanded  mode  for  the  data  part  of  record 
(array  A),  NUM  is  an  integer  number  equal  to  the  number  of  words  remain- 
ing in  the  record,  and  A is  an  array  of  variable  length  containing  the 
actual  data  to  be  read  or  written. 

The  logical  records  in  the  file  are  structured  as  follows 
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Tape  header 
Matrix  A header 
Column  1 of  matrix  A 
Column  2 of  matrix  A 

• 

Last  column  of  matrix  A 
Matrix  A trailer 
Matrix  B header 
Column  1 of  matrix  B 
Column  2 of  matrix  B 

Last  column  of  matrix  B 
Matrix  B trailer 


Tape  trailer 


The  contents  of  the  header  and  trailer  records  are  as  follows 


Tape  header 

IC0L 

-10 

K0DE 

0 

NUM 

7 

Data 

(TNAME(I) ,1*1,6) ,TM0D 

Tape  trailer 

-20 

0 

1 

0 

Matrix  header 

- 1 

0 

9 

(MNAME(I),I*1 ,6) ,MM0D,IMAX,JMAX 

Matrix  trailer 

- 2 

0 

1 

0 

where  TNAME  and  MNAME  are  six  character  tape  and  matrix  names  stored 
one  character  per  word  with  blank  word  fill  In  trailing  words  as  necess- 
ary, TM0O  and  W0O  are  Integer  numbers  used  as  modifiers  of  the  tape 
and  matrix  names  for  Identification,  and  IMAX  and  JMAX  are  the  row  and 
column  dimensions  of  the  matrix. 

In  matrix  column  data, IC0L  Is  a positive  Integer  equal  to  the  column 
number.  If  a column  of  a matrix  has  more  than  50X  non-zero  elements,  the 
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column  Is  written  In  expanded  form  (K0DE*O)  and  NUM  Is  equal  to  IMAX. 

If  a column  has  less  than  50X  non-zero  elements,  the  column  Is  written 
In  compressed  form  (K0DE»1).  Compression  Is  accomplished  by  forming  an 
array  of  the  non-zero  values  and  their  corresponding  row  locations  In 
the  column  as  the  data  part  of  the  record.  For  example.  If  the  12th 
column  of  a matrix  with  a row  dimension  of  100  had  only  3 non-zero  elements. 
It  would  be  compressed  and  appear  as  a logical  record  of  the  form 
IC0L  K0DE  NUM  Data 

12  1 6 Vl » Ll » V2’  L2’  V3’  L3 

where  IC0L  3 12  Is  the  column  number,  K0DE  3 1 Indicates  a data  record 
In  compressed  form,  NUM  3 6 is  the  number  of  words  remaining  in  the 
record  (three  pairs  of  value  and  location),  the  element  values  are 
Vg,  Vj,  and  the  row  locations  of  these  elements  are  , Lj,  and  L^. 

Note  that  1^,  L2,  L3  must  be  In  ascending  order. 

If  a column  of  a matrix  Is  completely  null,  no  record  for  that  column 
will  be  present  In  the  file. 

File  Contents 

The  files  shown  In  Figure  A.l  are  the  master  files  used  by  the  IMPACT 
system  to  pass  data  between  each  of  the  stand  alone  programs.  The  cont- 
ent of  each  of  the  card  Image  files  is  given  in  detail  in  Part  2.  The 
content  of  each  of  the  matrix  data  files  is  given  In  Table  A.l.  The 
FORMAT  tape  and  matrix  names  shown  In  parentheses  are  specified  by  user 
Input  and,  consequently,  may  vary  while  all  other  names  are  fixed  since 
they  are  Imbedded  In  the  program  code. 

USE  OF  UPOATE 

The  design  of  the  IMPACT  system  provides  for  the  optional  use  of  the 
COC  utility  UPOATE  during  the  preprocessing  phase.  This  phase  Includes 
the  execution  of  the  Laminate  Generator,  Initial  Generator,  and  Loads 
Generator.  Figure  A. 2 shows  the  Intended  use  of  UPOATE  In  an  application 
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TABLE  A.1 . CONTENT  OF  MASTER  MATRIX  DATA  FILES 


FILE 

10 

FORMAT  NAMES 

AND  MODIFIERS 

DESCRIPTION 

UZER0.1 

U . original  joint  coordinates 

MPT.l 

MPT,  material  property  data 

ECT.l 

ECT,  element  constants 

PRUPT.l 

PRUPT*  reorder',n9  transform  for 
degrees  of  freedom 

PRUF.l 

PDne»  reordering  transform  for 

RUf  elements 

KEL,1 

k,  lumped  element  stiffness 

MEL.l 

m,  lumped  element  mass 

Ml 

MTAPE.l 

KBEL.l 

k,  unlumped  element  stiffness 

CBEL.l 

c,  unlumped  element  damping 

FFBAR.l 

Fp,  element  force  transform 

SI6FB.1 

op,  element  stress  transform 

EPSIG.l 

c , element  strain  transform 

0 l 

DEBT.l 

<siy,  element  thermal  deformations 

EYT.l 

EYT,  element  variables 

C»NST,1 

C0NST,  problem  constants 

B 

LTAPE.l 

DPDT.l 

<P^j,  incremental  applied  loads 

I 


TABLE  A.l.  CONTENT  OF  MASTER  MATRIX  DATA  FILES  (Continued) 


FILE 

ID 

FORMAT  NAMES 

AND  MODIFIERS  • 

DESCRIPTION 

M3 

(FILE20.1) 

(PBUF.l) 

MPT.l  ^ 

UZER0 , 1 

ECT.l 

MEL.l 

KBEL.l 

CBEL.l  V 

FFBAR.l 

SIGFB.l 

EPSIG.l 

DEBT,! 

EVT.l 
CONST,! J 

Puc,  element  force  modal 
transformation 

Copied  from  file  Ml 

M4 

(FILE21 ,1 ) 

(PBPHI,! ) 

(PBUPTJ,! ) 

$P# transformed  Incremen- 
tal  applied  loads 

m 

P.ipT  i • degree  of  freedom  modal 

u J transformation 

< 

TABLE  A.1.  CONTENT  OF  MASTER  MATRIX  OATA  FILES  (Continued) 


FILE 

10 

FORMAT  NAMES 

ANO  MODIFIERS 

DESCRIPTION 

TAPE 

MATRIX 

C0NST.1 

C0NST,  augmented  problem 
constants 

TIME.l 

t,  Incremental  time  history 

) 

WIND0W, 52877 
(linear) 

UZER0.1 

(PBPHI.l) 

Copied  from  file  M3 

Copied  from  file  M4 

M5 

(PBUPTJ.l) 

Copied  from  file  M4 

TAPE.l 

(non-linear) 

RESPNS.8 

• 

V 5V  "8’  and  V 

i 

modal  response 

BARS, 6 

F„  , o.  , and  efl  , bar 
bB  *B  bB 

element  forces,  stresses, 
and  strains 

MEMBRN.B 

F^  * • an<*  ea  » membrane 

3|y|  ™ 

element  forces,  stresses, 
and  strains 

CELLS, 8 

F«  , dj  , and  c.  , cell 
^SC  SC  8C 

element  forces,  stresses, 
and  strains 

where  8 Is  the  Increment  number 
and  these  four  matrices  are 
repeated  for  each  time  Increment. 
If  there  are  no  elements  of  a 
given  type,  the  corresponding 
matrix  will  not  be  present. 
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using  all  three  preprocessing  programs.  The  basic  function  of  UPDATE 
In  the  IMPACT  system  Is  to  edit  the  model  definition  data  as  necessary 
for  each  of  the  preprocessing  programs.  A secondary  function  Is  to 
permanently  store  the  model  definition  data  In  an  easily  accessible 
form  once  It  has  been  read  Into  the  system.  On  a large  model  with  a 
large  amount  of  data,  this  will  avoid  many  card  handling  and  reader  errors 
each  time  the  model  definition  data  Is  required. 

In  Figure  A. 2,  each  of  UPDATE  steps  as  well  as  the  execution  of  the 
Laminate  and  Loads  Generators  are  optional.  Obviously,  many  combinations 
of  these  programs  are  possible  In  applying  these  tools  to  the  variety 
of  design  problems  a user  might  encounter.  Suffice  to  say  that  the  eff- 
ective use  of  UPOATE  In  the  preprocessing  phase  will  enhance  the  users 
effectiveness  by  alleviating  card  handling  problems. 
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SECTION  III 

OPERATIONAL  CONSIDERATIONS 


n 


GENERAL 

Only  one  of  the  six  new  programs  In  the  IMPACT  system  currently 
uses  the  overlay  capability  of  the  CDC  operating  system.  This  feature 
Is  used  extensively  by  the  Initial  Generator.  The  overlay  Is  accompl- 
ished using  the  SEGLOAD  option  of  the  loader.  The  SEGLOAD  directives 
for  Implementation  of  this  program  are  given  In  Section  V under  the 
heading  CORE  UTILIZATION. 

Significant  savings  In  core  requirements  could  be  realized  by  using 
overlay  In  the  Implementation  of  the  Linear  and  Non-Linear  Incremental 
Solution  programs.  The  remaining  new  programs  coded  for  IMPACT  are 
relatively  small  and  would  not  yield  enough  savings  to  justify  overlay. 
The  FORMAT  program  already  uses  this  feature. 

The  following  tables  gives  1)  minimum  core  requirements  for  the 
current  version  of  each  program  as  implemented,  and  2)  the  associated 
size  of  blank  common  In  the  compiled  code  as  Implemented  where  appli- 
cable. The  core  requirement  Is  In  octal  and  the  blank  common  region  Is 
In  decimal. 


PROGRAM 

CORE 

REQUIREMENT 

(octal) 

BLANK  COMMON 
LENGTH 
(decimal ) 

Laminate  Generator 

130000 

NA 

Initial  Generator 

130000 

NA 

Loads  Generator 

71000 

NA 

FORMAT  Phase  II 

100000 

10000 

Linear  Incremental  Solution 

104000 

8000 

Non-Linear  Incremental  Solution 

134000 

8000 

Postprocessor 

60000 

10000 
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EXTENDING  BLANK  COWON 


r w 


i 

i 


■ 


The  FORMAT  Phase  II  program,  the  Linear  and  Non-Linear  Incremental 
Solutions,  and  the  Postprocessor  use  blank  common  exclusively  for  core 
storage  of  data  during  their  execution.  Space  Is  allocated  within  the 
blank  connon  region  for  all  arrays  required  In  each  program.  The  array 
space  Is  allocated  according  to  problem  size  as  defined  by  Input  data. 

Even  though  the  blank  connon  regions  are  of  fixed  dimension  In  the 
source  code,  the  effective  blank  conrnon  region  can  be  Increased  at  exe- 
cution time  without  recompiling  any  code  when  operating  on  the  CDC 
system.  By  requesting  additional  core  on  the  job  request  card  and  supp- 
lying Input  data  to  the  program  Indicating  the  Increased  size  of  blank 
connon,  larger  problems  can  be  processed  than  could  be  accommodated  by 
the  dimensioned  size  of  the  work  area  In  the  compiled  code. 

For  example,  a program  requires  a field  length  of  100000  words 
(octal)  with  a blank  connon  region  of  10000  words  (decimal).  For  a 
given  problem,  an  Increase  of  2000  words  (decimal)  Is  required  In  the 
blank  connon  region.  By  requesting  an  additional' 4000  words  (octal) 
for  program  execution  to  account  for  the  2000  word  (decimal)  Increase 
In  blank  connon  region  and  Inputlng  data  to  the  program  specifying  this 
Increase  In  available  work  space,  the  program  could  process  the  larger 
problem  without  recompiling  any  code. 

Of  course,  this  procedure  can  be  used  only  with  a program  whose  str- 
ucture and  Input  are  consistent  with  this  philosophy.  The  four  programs 
of  the  IMPACT  system  mentioned  above  conform  to  these  requirements.  The 
blank  common  size  requirements  and  associated  Input  data  to  over  ride 
the  size  parameters  In  the  code  are  given  in  Sections  VII,  VIII,  and  IX 
of  this  document  under  the  heading  LIMITATIONS.  Again,  input  card  for- 
mats associated  with  this  data  are  given  In  Part  2.  Similar  Information 
regarding  FORMAT  Phase  II  Is  given  In  References  6 and  10  and  discussed 
In  Part  2. 
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OVER  RIDING  FILE  NAMES 


On  the  COC  system  when  executing  multiple  stand  alone  programs  as  may 
be  done  using  IMPACT,  It  may  be  necessary  to  over  ride  system  default 
file  names  because  of  conflicts  between  two  programs.  That  is,  one 
program  may  output  its  data  on  a file  named  TAPE!  and  the  second  program 
may  expect  Its  input  on  TAPE3.  In  order  to  resolve  this  conflict,  the 
default  file  name  In  one  of  the  programs  must  be  temporarily  renamed 
during  the  execution  of  that  program. 

On  the  CDC  system,  file  names  of  FORTRAN  programs  are  established 
by  means  of  the  PROGRAM  statement,  the  first  statement  of  the  main  pro- 
gram. The  system  provides  for  over  riding  these  names  at  execution  time 
by  means  of  a system  control  card.  The  format  of  this  command,  however, 
requires  foreknowledge  of  the  sequence  of  appearance  of  file  names  In 
the  PROGRAM  statement  for  a given  program. 

Table  A. 2 gives  the  PROGRAM  statements  for  each  of  the  seven  programs 
In  the  IMPACT  system.  It  shows  all  files  that  are  used  by  each  program, 
and,  together  with  Figure  A.l,  Impllclty  defines  all  files  used  as 
scratch  as  well  as  those  used  for  master  Input/output. 

SYSTEM  CONTROL  CARDS 

As  previously  stated,  the  seven  programs  of  the  IMPACT  system  may 
be  used  In  any  number  of  combinations  In  actual  application  to  the  var- 
iety of  design  problems  that  may  be  encountered.  However,- In  order  to 
show  an  example  of  a complete  deck  set  up  Including  system  control  cards 
and  data,  a single  run  using  all  components  of  the  IMPACT  system  Is  des- 
cribed here. 


Table  A. 3 gives  the  system  control  cards  for  a complete  solution 
while  Table  A. 4 lists  the  corresponding  user  Input  data.  The  complete 
job  Is  broken  down  Into  steps  for  reference  In  the  tables  and  this  dls- 
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TABLE  A. 2.  PROGRAM  AND  FILE  DECLARATIONS 


LAMINATE  GENERATOR 

PR0GRAM  LAMGEN  ( TAPE1 .TAPE2, OUTPUT, TAPE6«pUTPUT  ) 

INITIAL  GENERATOR 

PROGRAM  BIRDG1.  ( TAPEl -51 2, INPUT, OUTPUT, TAPE 5* INPUT, TAPE6-0UTPUT, 
TAPE?  -512,  TAPE8  -512,  TAPE9  -512,  TAPElO-512, 

TAPE11-512,  TAPE! 2-51 2,  TAPE! 3-51 2,  TAPEl 4-51 2, 

TAPE15-512,  TAPE! 6-51 2,  TAPE17-512,  TAPE! 8-51 2 ) 

LOADS  GENERATOR 

PROGRAM  L0DGEN  ( TAPEl .TAPE2 .TAPE3 .OUTPUT .TAPE6-0UTPUT  ) 

FORMAT 

PROGRAM  FORMAT  ( TAPE1-512,TAPE2«512,TAPE3-512,TAPE4-512, 

1 INPUT-51 2.TAPE5- INPUT, 0UTPUT-51 2 ,TAPE6-0UTPUT,TAPE7*51 2 , 

6 TAPE8-51 2 .TAPE9-51 2 .TAPE! 0*51 2 .TAPE! 4*51 2 .TAPEl 5-51 2 .TAPE! 6-51 2 ) 

LINEAR  INCREMENTAL  SOLUTION 

PROGRAM  RESPNS( TAPEl  -512.TAPE2  -512.TAPE3  -512.TAPE4  - 512 
1, INPUT  -512.TAPE5  * INPUT .OUTPUT *512, TAPE6-0UTPUT , TAPE7  -512 

2. TAPE8  -512.TAPE9  -51 2, TAPEl 0-51 2. TAPE11 -51 2, TAPEl 2-51 2 

3 . TAPEl 3-51 2 .TAPEl 4-51 2 .TAPEl 5-51 2 .TAPEl 6-51 2 .TAPEl 7-51 2 .TAPE! 8-51 2 

4 . TAPEl 9-51 2.TAPE20-51 2 .TAPE21-51 2 .TAPE22-51 2 .TAPE23-51 2 

5. TAPE29-512  ) 

NON-LINEAR  INCREMENTAL  SOLUTION 

PROGRAM  RESPNS (TAPEl  -512.TAPE2  -512.TAPE3  -512.TAPE4  -512 
1, INPUT  -512.TAPE5  -INPUT, OUTPUT- 51 2. TAPE6-0UTPUT, TAP E7  -512 

2. TAPE8  -512.TAPE9  -51 2. TAPEl 0-51 2, TAPE11 -51 2, TAPEl 2-51 2 

3 . TAPEl 3-51 2 .TAPEl 4-51 2 .TAPEl 5-51 2 .TAPEl 6-51 2 .TAPEl 7-51 2 .TAPEl 8-51 2 

4. TAPE1 9-51 2.TAPE20-51 2 .TAPE21-51 2 .TAPE22-51 2 .TAPE23-51 2 

5 . TAPE30-51 2 .TAPE31 -51 2 ) 

POSTPROCESSOR 

PROGRAM  POST  ( TAPEl .TAPE2.TAPE3, 

1 INPUT.TAPE5- INPUT, OUTPUT .TAPE6-0UXPUT  ) 
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TARLE  A. 3.  SYSTEM  CONTROL  CARDS  FOR  COMPLETE  SOLUTION 

STEP 

SYSTEM  CONTROL  CARDS 

m 

»C9,T07S,IOH<5,C*lU0  000#lTC3i.  D7SOS«5»*OM!I#  J20«fc 

Ll"ITC?7?7) 

i 

~UR0ATE(*,P,CpTAPE1) 

PETURN(NE«*»L> 

RE« I NO ( TAPED 

2 

ATTACH(IG0#LAmGENL60) 

*AP(ON) 

LCO. 

RETuRN(TapE1 .LOO) 

REwI^D ( T AP£2 ) 

3 

UPDATE(M,^,C»C*SE) 

RETURN(NE*PL,TAPf2) 

RE^INO(CARE) 

D 

ATTAO(B3t.CSD,8lRDG18L> 

UROATE(N.E) 

RC-IwOfCOHPILE) 

5 

ATTACH(8LG0,8!RD6as6) 

MAP(ON) 

SEGLQAO ( I*C0*,pILF#B«8aPS) 

L08eT(PRE3ET«ZERO) 

LOAD(BLGO) 

NOGO. 

RETuRm(B8LC»D#C0mPILE#BLG0,NEWpl) 

REWIMD(8AB8) 

8ABS(TA»E2,CASC) 

•ETu»M(8aR8) 

•EWIMP(T1PE2.CASE) 

6 

ATTACH? LGOtLOQGENL GO) 

PAP(QN) 

LCO(CASE,TX) 

■CTU«n(CA3E,Tx»LG0) 

P€wl*Df TAPES) 

7 

ATTACH(Ph2,phA5E2*,CY»1  ) 

»H2. 

•EwlN0(TAPf2,TAPpS,TAPF«) 

8 

•**2(,,tape«»,  taped 

RC*I*0(T  Ap£2#  TAP83#  TApf «) 

a 

“PH2(,,,TAPf20,,,,,TAPE21) 

•EWIM0JTAPE20.TAPE21) 

10 

ATTAC"(lG0*LlKlEA9LG0) 

HAP(QN) 

PPL ( 1 AAOOA ) 

100, 

RETURN (LGOt  TAPE 20*  T APf 21 ) 

_RfWlN0(TAP€2) 

11 

ATTAC*CLQO,ROSTI.GO) 

RAP  (ON) 

L80. 
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INPUT  DATA 


CCK 

cm 

T3.0 

5332 

1.0CM 

MM 

1 

3 

k.o 

5 

•.0 

7 

k.o 

k.o 

9 

• ,o 

11 

0,0 

k.o 

13 

12.0 

1« 

12.0 

15 

12.0 

k.o 

Ik 

12.0 

k.o 

MM 

1 

Ik 

MM 

5555 

1 

1 

1.0 

2 

1 

1.0 

3 

1 

a 

I 

1.0 

5 

2 

1.0 

k 

2 

7 

3 

1.0 

6 

3 

1.0 

5 

3 

10 

a 

1.0 

11 

a 

1.0 

12 

A 

13 

to 

2 § 

tk 

5555 

too 

101 

too  1 

1 1.0 

155 

5555 

l 

11  1.0 

5555 

100 

202k- 

Tkl 

1P1 

1 10.k 

t«k 

102 

1 7,kk 

(M 

103 

1 3*000, 

10k 

1 kkkoO, 

105 

1 .ok»* 

10k 

1 .33 

107 

1 .12k 

CM 

too 

1 .255 

C-3 

105 

1 50. 

110 

1 l.o 

e-o3 

155 

5555 

1 

II 

,00k 

3000 

2 

3 

3 

• 

•1.0 

5 

11.® 

2.0 

5 

1 

1.0 

k 

2 

k.o 

1.0 

k 

S 

0.0 

1.0 

MM 


1 


TABLE  A. 4.  INPUT  DATA  FOR  COMPLETE  SOLUTION  (Continued) 


INPUT  DATA 

il  i 

A 

11 

13 

ll 

tl  2 

I 

? 

A 

11 

11  3 

1 

3 

1 

T 

llffffffff 

10  t 

13 

1A 

too 

1.0 

20  2 

11 

14 

too 

1.0 

ioaaaa 

10  t 

A 

11 

11 

13 

too 

.1 

10  2 

10 

12 

14 

1A 

100 

.1 

10MM 

AO  1 

l 

3 

T 

1 

AO  2 

AOAAAA 

1 

T 

11 

A 

10  | 

A 

;* 

(•2 

.« 

e«i  .a 

£•» 

10  2 

1 l 

.« 

(•! 

.a 

£•1  .• 

£-1 

10  S 

13 

£•2 

.A 

£-2  .« 

£•2 

10  A 

AOAAAA 

11 

£•2 

.« 

e>2  .■ 

£•1 

•OCC*  4* SCI 
»tno 


Iff  3**  AT  JtanSa»6 

Iff  Ufa  CO,  LOGIC 

i**ut  t*»e  cAr**e,n 
iNffvit  ta*e  (LT**e,n 
OUTffUT  TA*E  <OecO"ff|t> 
IlHSTffUCTtON 

•T.OffT  a OffOT  .OEJCOL,  (1) 

•*vc  (orcoMff)  oaoT 

Mff  a 

ffffu*J,*auffS  a ffauffT  .oejCPL.  (««1 
*ff  a affuffj  ,"UlT.  at 

ffffT  a aff  .TffANga, 

*#LTL«U»  a •»U*»«eL  ffffT 

UT  a ffffuffT  ,T«ULT,  .jff 

SAVE  (OeCO«ff)  M*,AffUffJ.LTL.Uff 
•ff I*aT  (,,,J  UT 

1ENO 

1ffO*«AT  JTANQAaO 

JffUN  Gn,  LOGIC 

IlffuT  TiffJ  (HTAffCtt) 

txffuT  t»»e  (occo*’ff.n 

OUTffUT  TAfft  (T*A*«ff,l) 
IlNJTauCTlOh 

VAL, ra  a «a  .ulCffOfc,  ltl 

ffffTNT  {,,,1  VAL 
ffffuffj  a T*  ,T“ULT.  ffffuaj 

••MJff  a Tff  ,T«ULT,  ffffiia 

GffffO  a fftuffj  ,"ULT.  OAOT 
SAVE  (Tff AWJff ) TB,ffEUffJ.ff»'.iff,OMO 
•OOEI  • ff*UffJ  ,t»ult,  T* 

•• l N T 01)  "OOEI 

IlffECIAL 

» t«or«* 

UNO 


TABIE  A. 4.  INPUT  DATA  FOR  COMPLETE  SOLUTION  (Continued) 


sgognat 

ST  ANOANO 

GO,  LOGIC 

t*mjr  t*ge  (»ta*e,i) 

ISGtlT  T»(*e  (TGA*|0,l) 

outgut  r**c  oileeo, 

n 

OUTPUT  TAPE  (GlLlEl, 

n 

*1*»T»uCT!On 

UV( 

(GILE101  GSllG 

SAVE 

(EILEJ01  ngt 

SAVE 

(EtLEJOl  U*E*0 

Save 

(»!LEJO)  ECT.“EL»K8BL*CSEL.EEBA» 

SAVE 

(ETLEJO)  STGM, 

EGS16.0E8T.EVT, CONST 

SAVE 

f*ILEJl)  0090 

SAVE 

criietl)  0»UOJ 

SEnO 

1 

A tO  H» 

IG  18000 

10.0 

>0.0 

lo'n  ao;o 

80.0 

GO  , 0 

70. 0 

10.0 

00.0  100.0 

~Y~ 

1 

til 

1 1 

lit 

8 

lit 

G 

ttt 

7 

C AGO  t. 

GT»8T  TITLE  LINE 

• 

CAPO  I, 

GIG'ST  TITLE  LINE 

• 

CAPO  t. 

SECOnO  title  line 

to 

CAPO  I. 

secoNO  title  line 

>0 

GPSP 

1*  to 

JO 

t-  to 

cussion.  The  following  files  are  assumed  to  be  permanently  stored  disc 
files  In  the  system. 


LAMGENLG0 

BIRDG1SL 

BIRDG1LG0 
L0DGENLG0 
PHASE2A,  CY-1 
LINEARLG0 
P0STLG0 


Laminate  Generator  relocatable  decks 

Initial  Generator  SEGL0AD  directives  preceded 
by  a *DECK  UPOATE  control  card  (card  Images) 

Initial  Generator  relocatable  decks 

Load  Generator  relocatable  decks 

F0RMAT  Phase  II  absolute  file 

Linear  Incremental  Solution  relocatable  decks 

Postprocessor  relocatable  decks 


During  this  discussion.  Figure  A.l  should  be  referred  to  for  master 
input/output  files  of  each  program.  Also,  Table  A. 2 should  be  referred 
to  for  the  sequence  of  appearance  of  file  names  In  the  PR0GRAM  declara- 
tions of  each  program.  It  Is  assumed  that  the  reader  has  a working 
knowledge  of  the  CDC  system  and  Its  control  cards. 

In  Table  A. 3,  the  job  card  requests  140000  words  of  central  memory 
for  the  run  which  Is  more  than  sufficient  to  load  and  execute  all  prog- 
rams. The  LIMIT  card  requests  more  than  the  system  default  disc  storage 
space  normally  allocated  to  a job. 

The  CDC  utility  UPDATE  Is  executed  first  (Step  1)  receiving  input 
from  the  system  Input  file  and  outputlng  the  C0MPILE  file  under  the 
name  TAPE!.  In  Step  2,  the  Laminate  Generator  uses  all  default  file 
names  accepting  file  TAPE1  as  Input  and  creating  file  TAPE2  as  output. 
TAPE1  Is  the  partial  model  definition  data  and  TAPE2  Is  the  augmented 
model  definition  data.  This  Is  followed  by  another  execution  of  UPDATE 
(Step  3)  simulating  a final  editing  and/or  save  of  this  data  which,  rather 
than  being  output  under  the  default  file  name  C0MPILE,  Is  renamed  CASE. 

Step  4 executes  UPDATE  again  to  create  a card  Image  Input  stream 
for  SEGL0O  which  Is  output  under  the  default  name  C0MPILE. 


Step  5 Includes  both  the  loading  of  the  Initial  Generator  and  Its 
execution.  Input  directives  to  SEGL0AD  Is  the  card  Image  C0MPILE  file 
from  Step  4.  Input  to  the  Initial  Generator  Is  the  complete  model  defi- 
nition data,  card  Image  file  CASE.  The  normal  output  file  TAPE1  Is  over 
ridden  and  renamed  TAPE2. 

Next  the  Load  Generator  Is  executed  (Step  6).  The  normal  Input  file 
TAPE!  Is  renamed  CASE  In  order  to  accept  the  same  model  definition  data 
used  by  the  Initial  Generator.  The  output  Is  on  TAPE3,  the  default  file 
name  for  the  binary  matrix  data  output  file.  However,  the  normal  card 
Image  output  file  TAPE2  Is  renamed  to  the  dummy  name  TX  since  TAPE2  Is 
already  In  use  as  the  Initial  Generator  output  file. 

The  following  three  steps,  7,  8,  and  9,  execute  the  F0RMAT  Phase  II 
program  to  assemble  the  input  files  necessary  for  the  Linear  Incremental 
Solution  which  follows.  In  the  first  execution,  all  default  file  names 
are  used.  Input  Is  on  file  TAPE2  from  Step  5 and  TAPE3  from  Step  6. 
Output  Is  on  TAPE4,  In  the  second  execution  (Step  8),  Input  consists 
of  file  TAPE2,  the  same  file  Input  to  Step  7,  and  TAPE4  output  from  Step 
7.  In  this  step,  therefore,  the  normal  Input  file  name  TAPE3  Is  over 
ridden  and  renamed  TAPE4,  and  the  normal  output  fll^  name  TAPE4  Is  over 
ridden  and  renamed  TAPE3.  This  sets  up  the  stage  for  the  final  execution 
(Step  9)  In  which  the  Input  Is,  again,  file  TAPE2  from  Step  5,  and  file 
TAPE3  output  from  Step  8.  The  normal  output  file  names,  TAPE4  and  TAPE7, 
are  renamed  TAPE20  and  TAPE21  for  compatibility  with  the  following  step. 

The  Linear  Incremental  Solution  Is  then  executed  (Step  10)  using  all 
default  file  names.  Output  Is  on  TAPE2.  In  this  step,  blank  conmon  Is 
extended  to  accomnodate  the  size  of  the  problem.  On  the  first  Input  card 
for  this  step  (Table  3.3,  Step  10),  the  available  blank  cannon  Is  speci- 
fied to  the  program  In  decimal  as  15000  (last  value  on  card).  This  Is 


execution,  an  RFL  command  Is  placed  before  the  load  and  go  command, 

LG0.  This  effectively  over  rides  the  field  length  the  loader  normally 
establishes  from  the  relocatable  file  It  accesses  for  the  load.  The 
requested  field  length  on  the  RFL  card,  140000  words  (octal).  Is  suffi- 
cient to  compensate  for  the  7000  words  (decimal)  added  to  blank  common. 
(See  minimum  core  requirements  and  dimensioned  blank  common  block  sizes 
discussed  at  the  beginning  of  this  section). 


The  last  step,  execution  of  the  Postprocessor  (Step  11)  Is  then  per- 
formed using  all  default  file  names,  in  which  case,  TAPE2  is  the  input 
file. 


This  example  would  be  Identical  if  the  Non-Linear  Incremental  Solu- 
tion had  been  used  ra’ther  than  the  linear  with  the  following  exception. 
The  default  output  file  for  the  non-linear  Is  TAPE23.  The  last  control 
card,  the  command  executing  the  Postprocessor,  would,  therefore,  have 
to  over  ride  the  default  Input  file  name  of  TAPE2  to  TAPE23  and  would 
appear  as 

LG0  (,  TAPE23) 
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